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Abstract 

We match the Higgs sector of the most general flavour breaking and CP violating minimal su- 
persymmetric standard model (MSSM) onto a generic two-Higgs-doublet model, paying special 
attention to the definition of tan P in the effective theory. In particular no tan /3-enhanced loop 
corrections appear in the relation to tan f3 defined in the DR scheme in the MSSM. The cor- 
rections to the Higgs-mediated flavour-changing amplitudes which result from this matching are 
especially relevant for the Ba and Bg mass differences AMd^s for minimal flavour violation, 
where the superficially leading contribution vanishes. We give a symmetry argument to explain 
this cancellation and perform a systematic study of all Higgs-mediated effects, including Higgs 
loops. The corrections to AMg are at most 7% for /i > and Ma < 600 GeV if constraints 
from other observables are taken into account. For fi < they can be larger, but are always less 
than about 20%. Contrary to recent claims we do not find numerically large contributions here, 
nor do we find any tan /3-enhanced contributions from loop corrections to the Higgs potential in 
B^ r+z/ or 5 — > Xg'j. We further update supersymmetric loop corrections to the Yukawa 
couplings, where we include all possible CP-violating phases and correct errors in the literature. 
The possible presence of CP-violating phases generated by Higgs exchange diagrams is briefly 
discussed as well. Finally we provide improved values for the bag factors P^^^, P2^, and P^ll 
at the electroweak scale. 
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1 Introduction 

Supersymmetry constrains the structure of the Yukawa coupUngs of the minimal supersymmet- 
ric standard model (MSSM) to those of a special two-Higgs-doublet model (2HDM). In this 
2HDM of type II one Higgs doublet, only couples to up-type fermions, while the other one, 
Hd, only couples to down-type fermions. As a consequence, there are no dangerous tree-level 
flavour-changing neutral current (FCNC) couplings of the neutral Higgs bosons. However, the 
presence of supersymmetry-breaking terms destroys this pattern at the one-loop level, permitting 
couplings of both Higgs doublets to all fermions. Thus the resulting Higgs sector is that of a 
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Figure 1: Leading contributions to B^j — Bg mixing from supersymmetric Higgs bosons. The 
FCNC couplings are induced by supersymmetric loops. The coefficient of Qf^^ = {hRqL)(pR(lL) 
vanishes, if the tree-level relations between Higgs masses and mixing angles are used. 



general 2HDM, often called 2HDM of type III. As pointed out first by Hall, Rattazzi and Sarid, 
the loop-induced Yukawa couplings can compete with the tree-level ones in the limit of a large 
tan/3 = Vu/vd, which is the ratio of the vacuum expectation values (vevs) of Hu and Hd [1]: 
in the relationship between if„ ^-couplings and observed masses of the down-type fermions the 
loop suppression factor ~ 0.02 is offset by a factor of tan/?, so that 0{\) corrections to the 
type-II 2HDM are possible for tan/3 ~ 50. In such scenarios also loop-induced FCNC 
couplings of neutral Higgs bosons appear [2], which allow the branching fractions of (yet un- 
observed) leptonic B decays to exceed their standard-model values by more than two orders of 
magnitude [3]. This observation has stimulated a large activity in flavour physics and powerful 
constraints on the MSSM Higgs sector in scenarios with large tan [3 have been derived from 
B factory data [3-6]. These Higgs-induced effects in flavour physics are very transparent in 
the limit 

MsusY > Ma ~ , (1) 

where Msusy denotes the generic mass scale of the superpartners and the masses Ma, Mh, 
and Mh± of the five physical Higgs bosons are taken to be of the order of the electroweak scale 
V = ^/vl + v"^ = 246 GeV. All low-energy observables can be computed in the type-Ill 2HDM, 
which emerges as the effective theory in the limit of Eq. (1). The new couplings can be calculated 
from finite one-loop diagrams with supersymmetric particles and thus become functions of the 
MSSM parameters, so that the desired constraints on the supersymmetric parameter space can be 
derived. The effective 2HDM Lagrangian efficiently incorporates all large-tan P effects, equiv- 
alent to a perturbative all-order resummation of those radiative corrections which are enhanced 
by a factor of tan /3 [7]. 

Bg — Bg mixing (with q = d or s) plays a special role among the FCNC transitions of B 
mesons. Here the leading new effect stems from effective tree-level diagrams with neutral Higgs 
bosons (see Fig. 1). A priori the dominant contribution is expected from Yukawa couplings to 
right-handed b quarks, generating the effective AB = 2 operator 

Qf L ^ _ (2) 

However, the corresponding coefficient Cf^^ vanishes exactly, if one employs the tree-level rela- 
tions between the Higgs masses and mixing angles [2]. Nevertheless, sizeable effects in Bg—Bg 
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mixing are possible even in scenarios with minimal flavour violation (MFV) [8-16], in which 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [17] is the only source of flavour violation: 
keeping the strange Yukawa coupling non-zero one finds a non-vanishing contribution to the 
coefficient of 

= (hmL) {hqR} , (3) 

which depletes the Bg—Bg mass difference AMg [5]. The tree-level vanishing of Cf^^ calls for a 
systematic analysis of all subleading effects. In particular, the contribution that stems from Qf^^ 
can a priori compete with the contribution of the operator (^2^ above if the one-loop corrections 
to the MSSM Higgs potential [18-24] are taken into account. While a lot of work has been 
devoted to the analysis of the Yukawa sector [2, 3, 5-7], little attention has been given to effects 
from the Higgs potential. An exception is Ref. [25], which finds large contributions. We revisit 
these effects in the present paper and perform a systematic matching of the MSSM Higgs sector 
onto the type-III 2HDM. The result is not only relevant for the calculation of Cf^^, it also clarifies 
the relationship between the definitions of tan (3 in the MSSM and the effective 2HDM. This is 
important to link the constraints from flavour physics to other fields of MSSM phenomenology, 
in particular Higgs physics. Our paper is organized as follows. We derive the corrected B — B 
mixing amplitude in Sect. 2, including all relevant subleading contributions. The renormalization 
of tan (3 and som further technial issues are the subject of Sect. 3. In Sect. 4 we apply our new 
formulae to the phenomenology of B — B mixing, analysing the mass differences AMd and 
AMg as well as CP- violation. Our results are summarized in Sect. 5. We list our notation and 
our technical results in four appendices. Parts of our results were previously presented by one of 
us at a conference [38]. 



2 Higgs-mediated effects in B — B mixing 

The quantity governing the Bq — Bq mass difference is the off-diagonal element of the Bq — Bq 
meson mass matrix: AMq — 2 |M|]^|, with 

^ {B,\Ki=''\Bq) 

2Mb, 

The AB — 2 effective weak Hamiltonian H^^'^ consists in general of eight dimension-six 
operators: 

K^/^' - ^^^^ E C^{^^^)Q^i^^n) , (5) 

i=l 

with Xqf, = VtqVfl. The set of operators in Eq. (5) comprises the standard-model operator, 

Q^'-'' = (pLl.qL) (^Ll^qi) , (6) 
the two scalar operators defined in Eqs. (2) and (3), the operator 

gf ^ = {hqu) {hqu) , (7) 
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and four other operators. The complete list of operators plus the relevant evanescent operators 
is given in Eq. (135) and Eq. (136) of Appendix C. We express our results in terms of ma- 
trix elements at the high scale which we choose equal to the top mass fntijfit) = 164 GeV. 
In this way the other four operators do not appear in our formulae. However, some of them 
are needed to connect Qi(A*/i) with Qi(A*6) at the low scale /i^ ~ at which their matrix el- 
ements are computed, because they mix with Q^^, or under renormalization. We 
follow the conventions of Refs. [26] and [27] for operators and matrix elements. In particular we 
parametrize the hadronic matrix elements as 

{B,\Q,{^u)\B,) = '^-MljlP,. (8) 

The Pj's are obtained [27] by renormalization-group evolution from the conventional bag factors 
Bi computed at the low scale ^5. We calculate the Pj's from up-to-date lattice QCD results in 
Appendix C, where we fully exploit constraints from heavy-quark relations. This is a new feature 
of our analysis compared to previous studies of new-physics effects mB — B mixing. 



2.1 Effective tree-level Higgs exchange 

The Higgs sector of the MSSM contains two SU(2) doublets and H^, 

ht\ fhiV ( 1 



u 



of hypercharge -1-1/2 and —1/2, respectively, with vacuum expectation values (vevs) (/i^ = 
^u,d/ of relative size tan /3 = Vuj v^. Integrating out supersymmetric particles, the Lagrangian 
of the resulting effective 2HDM is no longer restricted to be of type 11, and is constrained only 
by the electroweak symmetry. Neither will it be renormalizable, with operators of dimension 
greater than four encoding effects that decouple at least as w/Msusy for heavy superpartners. 
We begin with a short review of some pertinent aspects of the general 2HDM. 
Defining 

■ $ \ _ / cos/3 sin/3 \ / -eEl 
^' )~\ - sin/3 cos/3 ) \ 

the most general fermion-Higgs interactions up to dimension four read 



(10) 



Cy^ dm Md^^ ^^Qlj - dm Kij ^'^Qlj 

- —'^Ri ^Uii Qlj -^-URi KijQLj + h.C. , 

where we have employed the notation a ■ b = a^e b. By construction, the vev of $' vanishes, 
whereas $ has ($) = (0, v/v^)^ and contains all three Goldstone bosons. Hence only $ can 
contribute to the fermion masses and only can have flavour-violating neutral couplings. The 
flavour basis is defined such that the down-quark mass matrix is diagonal. In this basis the 
FCNC Higgs couplings to 6-quarks are governed by n^q or Hqi, (q — dor s). 
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The renormalizable Higgs self -interactions are comprised in the most general gauge-invariant 
dimension-four two-Higgs-doublet potential [28], 

V = ml^HlHd + ml^HlH^ + {ml^Hu ■ Ha + h.c.} 

+ ^{HlHaf + ^{HlH^f + \,{HlH^){HlHa) + X,{HlHa){HlH^) 

+ 1^ • Haf - HHlHa) {H^ ■ Ha) - HHIH^) {H^ ■ Ha) + h.c.j . (12) 

The couplings m\2, A5, Ae, and A7 are in general complex, yet the vevs v^.d can be made real 
by a [/(I) transformation on the Higgs fields. The definitions of mf^ and Aj in Eq. (12) coincide 
with Ref. [28] except for A3 and A4: we associate a different operator with A4 to eliminate it from 
tree-level neutral-Higgs phenomenology and have instead A3 = \^^^ + Xf'^^ and A4 = —X^^^. 
Shifting the fields in Eq. (12) by their vevs, which minimize V at tree level, ^ 

K,d = '^^'"^'d + 'i^^^d + iXu,d) , (13) 

determines the physical Higgs-boson mass matrices and interactions. We write the neutral-Higgs 
mass matrix in the basis (0^, Xd: Xu) in terms of 2 x 2 blocks, 

^ / Ml MIj 
Ml^\ 1 , (14) 

with M^, M'^j, and M| given in Eqs. (23-26) below. In the CP-conserving case, M|j = 0, and 
Af|. and M| are diagonalized by rotating the CP-even and CP-odd Higgs fields through angles a 
and 13, respectively: 

(t>d \ ^ ( cosa -sina \ f H^ \ f Xd \ ^ f cos/3 -sin/3 \ f G° \ ..a. 
(t)u ) V sina cosa ) \ hP ) ' \ Xu ) V sin/3 cos/3 ) \ ) ' ^ ^ 

The same angle /3 — oxcidiXiVu/va as defined above appears because (and only when) Vu, va 
minimize V. If CP violation is present, four physical mixing angles ai^2,z and /3 are required to 
diagonalize Mq . The charged-Higgs mass matrix is always diagonalized by /3, 

K\_( ^os(3 -sin/3\ /G+A 

K )~\ sin/3 cos/3 ; V ^+ ; ■ ^ ^ 

The non-standard effective operators Q\^, Qi^^, and Q'^-^ are generated at tree level via the 
exchange of neutral Higgs bosons (see Fig. 1) with the Wilson coefficients 



Tree level" here refers to the 2HDM. We defer a discussion of quantum corrections to w„ and V4 to Sect. 3. 
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and Cf^^ obtained from Cf^^ through the replacement {ntg f T~ {i^*qbY ^''*- We find that, 
in the general case, the Higgs propagation factors can be expressed as follows: 

_ det (M| + M| + zM|, - zMlJ) _ det B 

J Q<^Q — 999' \^^/ 

_ det {Ml - M] - iMl, - iMM) _ det A 

111 — 222' \^"/ 

where the denominators contain the product of the three nonzero eigenvalues of Mq. In the 
CP-conserving case, Eqs. (18) and (19) reduce to the well-known expressions 

± _ sin'(Q; - /3) cos^(a - ^) J_ 

Ml + Ml ^Ml' ^^"^ 

where Ma denotes the CP-odd Higgs-boson mass. 

The discussion so far has been completely general. Particularizing to the MSSM, a perturba- 
tive matching calculation relates the two theories. At tree level this trivially results in 

Mf = ^ cos /3 Fd, M(°) = ^ sin /3 F„, 

= - sin /3 Y^, k^^^ = cos /3 

mlT = l/^r + ml, ^ ml A?) = Af = -Af = + ,-)/4 ^ ^V4, ^^l) 

2(0) I |2 , 2 _ 2 \{0) 2 /o 

mlT = 5;., Af = Af = A?) = 0. 

At this order k^^^ and k^^^ are aligned with Mf^ and m1°\ respectively, so that no FCNC are in- 
duced, as it must be in a model II. At one loop, all couplings in Eq. (12) are generated. Moreover, 
the corrections to the Yukawa couplings have the more general form 

M^^^ = cos /3 [AYd + tan (3 AK] , k^^^ = - sin /3 [AYa - cot (3AK] , (22) 
v2 

where AY^^^ and AKij parametrize the one-loop vertices djuHa ■ Qlj and djuHlQij, respec- 
tively. Diagonalizing M^ rotates k^^\ giving rise to a flavour-violating coupling oc Y^ tan /3/(16 tt^), 
which can be of 0{1) for tan f3 ~ 50. 

The origin of this explicit tan (3 enhancement (in addition to the mere presence of large down- 
type Yukawa couplings), which can compensate the loop factor 1/(16 tt^), is the replacement 
of Vd by Vu ^ Vd in the contribution of AK to M^ [I]} This removal of a v^, suppression 
can happen only in dimensionful quantities. In the fermion mass terms, only one power of 



■^We tacitly assume that the fermion kinetic terms in the effective 2HDM have been made canonical. Such a 
field renormalization does not contribute factors of tan (3 because it is determined by dimensionless couplings. Cf . 
Sect. 3 for a discussion of field renormalization. Our AK and AY^ correspond to AuY^, and — A^Vd. respectively, 
in the first paper of Ref. [5]. 
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tan P can appear because there was only one power of Vd to begin with. This is in agreement 
with the findings in [7]. Our approach using un-shifted Higgs fields ("unbroken-theory") makes 

particularly evident that this result holds to all orders, as the Yukawa Lagrangian only involves 
dimensionless couplings and there are no hidden factors of tan (3. Although we have integrated 
out only the sparticles - as we assume a hierarchy v, Ma ^ Afsusv - the argument continues 
to hold if we also integrate over the Higgs fields, keeping only constant background values of 
(spurious). The reason is that for determining the mass matrices, the relevant external 
four- momenta are of 0{mg), providing an expansion parameter rriq/v or rriq/MA. Hence the 
Higgs contributions to the effective potential (which on general grounds respects the electroweak 
symmetry) can be organized into a (local) effective Lagrangian, with m^-suppressed corrections 
to the form Eq. (11) encoded in higher-dimensional operators with additional derivatives acting 
on dui or Qlj. The contribution from both Higgs and sparticle loops to M^ is then simply 
obtained upon substituting for $, $' their vacuum expectation values. This mass matrix is to be 
identified with a short-distance (such as MS) mass in the effective QCD x QED at low energies, 
where the dependence on the chosen scheme cancels against the explicit form of the matching 
(of the 2HDM onto QCD x QED). 

There is only one other place where a similar tan /3 enhancement can occur, namely in the 
dimensionful self-couplings of the (shifted) Higgs fields, that is, their masses and trilinear cou- 
plings. Indeed, at dimension two it is exhibited in the neutral Higgs mass matrix Eq. (14). 
Explicitly, one has (with sp = sin (3, c/s = cos (3, and = Re A^) 

m2 _ „2 f + '^y^s^cp + AicJ X\sl + X^spcp + A^c2 \ 2 

M] = Ml ( ± -'f^ ] , (24) 



where {m\2y has been traded for M\, with 



sp cp Ml = {ml,y - ^(A^, 4 + 2X1 spc^ + A^ c|). (25) 

If CP is conserved, in the limit of infinite tan (3 {cp — > 0) the leading mass splitting Mjj — M\ — 
As^;^, and the leading correction to the tree-level result a = is determined by Ay. In the former 
case, an enhancement by two powers of tan/3 occurs {Mfj — M\ = 0(cos^/9) at tree level), 
while the loop correction to a is enhanced by a single power of tan P with respect to its tree- 
level value. Either effect is sufficient to remove the cancellation in in Eq. (20). Moreover, a 
1 / tan /?-unsuppressed CP-violating contribution proportional to A| and X\ appears to occur: 

2 _vW Xi4 + 2Xlspcp -Xis^cp - 2Xic} \ 
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2 V 2A'7 sj + XI spc/s -2X\ S/scp - XI 

where Xf, = Im A^. However, as we show in Sect. 2.3 below, the individual phases of A5 and 
A7 become unphysical in the limit tan /3 — > 00, and mixing between the CP-even and CP-odd 
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sectors is described by a single angle ol , determined by the relative phase of A5 and Ay. Finally, 
the charged Higgs mass matrix is given by 



Here no tan enhancement due to loop-induced couplings occurs. 

Unlike the case of the fermion mass matrix, the typical momentum flowing through the ef- 
fective Lagrangian Eq. (12) for an on-shell Higgs is itself of 0{v) or 0{Ma). Hence Higgs-loop 
contributions to the Higgs masses cannot be included in Eq. (12), but rather the full effective 
action would be needed. Higgs-loop effects in and Kqb multiplying could, however, be 
included via Eq. (11), since again the momenta flowing through the vertices are much smaller 
than V, Ma- This is not possible in Higgs boxes, where large momenta flow through the FCNC 
vertices. We will present a systematic method to include all Higgs-loop contributions in Sect. 2.3. 

It is instructive to consider the explicit form of the numerator in Eq. (19), which is 



With Eq. (21), det A = ?/(AiA2 — A|) s^c^ = 0, reproducing the known vanishing of 
employing the tree-level MSSM Higgs sector. The cancellation is removed already at the leading- 
logarithmic level. For instance, A2 alone receives a large additive correction oc due to top- 
quark loops, which is also responsible for the most important correction to the tree-level mass 
of h. The corresponding corrections could be computed by RG-evolving the tree-level couplings 
in the effective 2HDM. However, as we are considering large tan /3, we expect (and find below) 
the most important effect to be due to A5 and A7, which remove the 0(c|) suppression of the 
leading-log result, as anticipated above. 

2.2 The case of minimal flavour violation 

From the discussion so far it follows that = 0{1/Ml) > \J^-\ = 0{l/{lQ7r^Ml)), 

implying |C^^| ^ jCf'^^l for generic k^,^ such that the motivation to consider JF^ at all is 
not very strong. The situation is fundamentally different for MFV because then the contribution 
proportional to turns out to be suppressed by a light quark mass, introducing a further small 
parameter mg/mft comparable to l/(167r^) or 1/ tan/3 for g = s (and negligible for g = d). For 
simplicity, in this paper we consider the simplest version of MFV, assuming flavour-universal 
soft breaking terms rfiQ, ml and mfi and trilinear SUSY-breaking terms T^.^., T^.^. which are 

^In Ref. [13], an argument based on SU{2) x f/(l) gauge invariance was used to infer that (in the present 
notation).?"" = 0{v'^ /M\). This statement, which clearly is respected by our Eq. (19) in conjunction with Eq. (28) 
(recall = 0{v'^)), is about the asymptotic behaviour as u/M^a — > 0. The latter is not necessarily a small number 
in practice. Indeed, many of the analyses in the literature have dealt with the case Ma ^ 200 GeV. 




(27) 



det^ 



v"" [{X2XI - xf) 4 + 2 (A2A; - A3A; + a;a7) 

+ (A1A2 - A^ + IA5I' - 2A6A; + 4A*A7) slcl 



+ 2(A5A; - A3A6 + XiX7)spd + (A1A5 - A^)ct 



(28) 
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proportional to the Yukawa matrices and therefore diagonal in the super-CKM basis (denoted 
with a hat): T^^. — atVui^ij and Td^. — aby^iSij, see Appendix A for details of our notation. 
The structure of our results, however, does not depend on these additional assumptions. The 
tan /5-enhanced loop-induced FCNC couplings of the neutral Higgs bosons in Eq. (11) can be 
expressed as: 



2 V2mb 
i^bq — Vt \b 



1 



1 



V cos^ (3 1 + €3 tan j31 + Sq tan j3 ' 



2 \ * 



V2 



rria 



V cos^ /9 1 + £3 tan /? 1 + eo tan j3 ' 



(29) 
(30) 



with ijt = V^rrit/ {v sin /3) and Xgb = VtqVtl- The effective couplings ey, cq and €3, which depend 
on the MSSM parameters, have been analysed in the decoupling limit Msusy ^ v in the limit 
g — g' = in Refs. [2-4] for the case that ey, eo and e-^ are real. We consider the general 
case allowing n, the universal trilinear term at and the gaugino mass parameters to be complex. 
Effects from non-zero g, g' have been taken into account in Ref. [5], where also effects beyond 
the decoupling limit were considered. The corresponding expressions for Msusy ^ v, suited 
for our analysis, were derived in Ref. [25]. We have recalculated the FCNC couplings of neutral 
Higgs bosons including all CP-violating phases and found agreement with the results for the 
FCNC self -energies given in Ref. [5], but encountered a significant discrepancy with Ref. [25]. 
In our results, the phase conventions of the first five parameters can be inferred from Eqs. (98) 
and (100) of Sect. A. The phase convention for M3 complies with that of Mi 2 and the gluino 
mass equals Mg = \Ms\. Of course one can choose one of these parameters (e.g. M3) real. Now 
the effective couplings of Eqs. (29) and (30) read: 

£0 = 



+ 



+ 



37r M3 



967r2 Ml 

9" ^ 
1447r2 Ml 



—1 a 

167r2 




+ 2H2 



327r2M2 




+ ^Y,v/M , 



£3 = eo + ey 



Here 



H2{x,y) = 



X logx 



ylogy 



X] 



(31) 

(32) 
(33) 

(34) 



[l-x){x-y) il-y)iy 

Numerically, the electroweak contributions in eo can be of 0(10%). They improve the compar- 
ison with the results computed with full chargino and squark mass matrices (see Eq. (5.1) in the 
second paper in [5]). 
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Ref. [5] also discusses threshold corrections to the fermion kinetic operators (wave function 
renormalizations). While these terms are not tan /3-enhanced, the flavour-diagonal quark wave 
function renormalization constants receive sizable contributions from squark-gluino loops. One 
can parametrize these loops in terms of a new quantity colkiu which will add to eg in the relation 
between the MSSM Yukawa coupling y^. and the physical quark mass m^. (see Eq. (103) for the 
case of the bottom Yukawa coupling). eo|kin will likewise appear in the relation between Kij and 
ydi, but it drops out once Kij is expressed in terms of md^, so that it does not appear in Eqs. (29) 
and (30). This cancellation of the flavour-diagonal quark wave function renormalization can be 
verified by inserting Eq. (2.29) into Eq. (2.26) of the second paper in [5]. This feature can be 
traced back to the fact that the wave function renormalization affects both the tree-level and the 
loop-induced Yukawa couplings with the same multiplicative factor. 

Comparing our result with Ref. [25], we find different results for eo and ey'. In Ref. [25] 
the chargino-stop contribution proportional to g"^ is erroneously assigned to ey rather than eq. 
Since this piece does not contain any Yukawa couplings (the chargino is a pure wino here), all 
three generations contribute in the same way and the resulting overall CKM structure combines 
to V*^yuq + V*^Vcq + V^^Vtq, wMch Is zero fox q ^ h and equal to one for q ^ b. This GIM 
cancellation eliminates the wino-stop loop from ey, while this loop contributes to eo twice as 
much as the corresponding loop with a neutral wino-like neutralino and a sbottom. The two 
terms are combined into the last term in Eq. (31). Omitting the chargino loop here would violate 
SU(2) gauge symmetry, which also enforces = 6^ in the decoupling limit. Since ey normalises 
all Higgs-induced FCNC couplings, one should verify the accuracy of the Msusy ^ v limit: It 
is easy to include the tan /3 -enhanced contributions to ey to all orders in i>/Msusy- To this 
end one merely has to calculate the FCNC bRqL self-energy using the exact chargino and up- 
squark mass eigenstates. This self-energy renormalises the off-diagonal pieces of the quark mass 
matrix and cause the mismatch between the flavour structures of the latter with the Yukawa 
couplings leading to ey ^ 0. In higher loop-orders tan /3-enhanced contributions are suppressed 
by products of small CKM elements (and are negligible) or are flavour-conserving and therefore 
contribute to eo rather than to ey . Using the buqi self-energy (S^^) (with q — di) from Ref. [5] 
one finds 

=J-<H, i'^, + ^ Mil! (35) 

(Note that (S^^^)'^* oc and be aware of the different sign conventions for ijij in Eq. (103) and 
Ref. [5].) We stress that Eq. (35) must be evaluated for i 7^ 3, so that the GIM cancellation of the 
above-mentioned wino-stop loop takes place. Numerically one finds a marginal impact of eY,v/M- 
Setting all supersymmetric massive parameters equal to a common value Msusy> one finds that 
ey,D/M amounts to a mere 1.4% correction to ey for Msusy = 400 GeV. Even for Msusy = 
150 GeV, for which the expansion in v/ Msusy formally breaks down, eY.v/M depletes ey by as 
little as 8%. eY,v/M also enters eo through Eq. (32). It can be inferred from Ref. [7] that this 
procedure indeed leads to the correct all-order resummation of the tan /3-enhanced corrections 
involving y^. Corrections to eo beyond the Msusy ^ v limit from g, g' and are considered in 
Refs. [7] and [5]. We remark that no terms proportional to yl occur in Eqs. (31-33), because the 
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corresponding loops violate hypercharge and involve a suppression factor of v^/Mgusv 

We verify from Eqs. (29) and (30) that n*f,Hbq multiplying ^+ in Eq. (17) is suppressed 
by a factor rriq/mij relative to k^^, which multiplies JF~. Hence Cf^^ is naively leading (over 
C^^) from the point of view of MFV alone, and a meaningful analysis of — Bg mixing 
requires a systematic investigation of all leading corrections to its vanishing "tree" value. (The 
coefficient Cf^^ both undergoes a strong my ml suppression anJ involves and can thus be 
disregarded.) It is then useful to think of the AF = 2 amplitude as being a function of the four 
small parameters identified so far: 

[A-kY mb tan/3 Msusy 

The vanishing 2HDM tree diagram for J^' is (superficially) 0{{cot p)''^ l'^ p'^ cu^), i.e. 0(1) 
when treating all expansion parameters on the same footing. Conversely, JF+ is nonzero at the 
tree level but is suppressed by one power of u, which is non-negligible only for q = s. We 
have already seen that JF" vanishes exactly for tree-level matching (or up to 0(1/ tan^ /3) when 
including leading logs), so there are no 0(1/ tan/3) corrections at first subleading order. This 
leaves loop corrections (via sparticle corrections to the Aj as well as loops in the effective 2HDM) 
and possible corrections due to higher-dimensional operators, not written in Eqs. (11) and (12). 
We now discuss these contributions in turn. 



Sparticle loops One-loop contributions from higgsinos, gauginos, and sfermions correct the 
values of Ai,2,3.4 in Eq. (12) and induce non-zero couplings Xq^j. As a technical result of our 
paper, we have computed the Aj for general sparticle masses and flavour structure. These results 
are reported in Appendix B. At tree-level in the effective theory and in the leading order of 
1/ tan/3 the quantitiy JF" receives only contributions from A2, A5, and A7, cf. Eq. (28). The 
general results of Eqs. (71),(116),(120),(121),(123), (125-126), and (128) for the MFV case read 

A7 = A7 ^Y^ji*""^'^^'^ (2g''' Co {mi, me, The) + (^fi^^ - Co{me,mi,mij^ 

+ fj,ar\^i\'^\yr\^Do {me,me,mi,mi) 
1 



+ (3fiM2D2 (IM2I, IM2I, H |/x|) - ^//|M2|S^ (IM2I, |//|)^ 



1 



^g'g'"fii^\M,\B', {\M,\, |//|) + (IM2I, |//|) 

-4(Mi + M2)L)2(|Mi|,|M2|,|/x|,|/x|) 

+ / (^fxM^D2 (|Mi|, |Mi|, H I/.I) - ^fx\MM i\Mi\, |/x|)) 



(37) 
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Afi — Ar 



16^2 



fi <3ai^\yb\ Do{md,md,mQ,mQ) + 3a^\yt\ Do {mQ,mQ,mu, tUu) 



+ 4 \yr\ ^Do (me, me, fhurhi) - WM^Do (IM2I , IM2I , |At|, 

- 2g^/M^M2Do {\M,\ , IM2I , - / MID^ {\M,\ , |Mi| , |/.|, |/.|) 



and 



(38) 



' ' 4 167r2 



- |5''*^o {rhe, rhe)-^(g^ + f ) 5o {mi,rhi) 



+ ^ (s-'" - l/uyr^Co {me,mi,mi) - g'''\nyr\'^Co {mi,me,rhe) 

- lnyrfDo {me,me,mi,mi) 

- ^9'* Bo {fhd.rhd) + (-3|yt|^ + '2g''\yt\^ - g'"^ Bo (m„,m„) 
+ ^ (-9^^ - 2A\y,\' - / - 4\yt\' (/ - 3/)) Bo {rriQ, mg) 

+ ^ (3/ - 12|yt|^ \atyt\'^Co{mQ,mQ,mu) 

+ 2 -3||/t|^j |at2/t|^Co (mQ,m„,m„) - g''\iJ.yb\'^Co {mQ,md,md) 

- ^ (3/ + /) l/^yfcl^Co (md, mg, ttiq) 

- ^atytl^Do {rfiQ, mg, m„, m„) - 3|//y6|^Do ("^d, rhq, fho) 

+ (31/^1/6^50 {^d^rhq) + lnyrl^B^ {rhe,rhi) + 3|at2/tp5o (m„, mg)) 

+ - -21og-|^' -61og-|(7 -81og-^^' 



/^o 



/^o 



/^o 



31og4(,^ + /)-log!^(9,^ + / 



~ 2 
7717 



1 2 /2 



24^ 



-12L'2 (IM2I, IM2I, |//|, |//|) |M2|2 - 60^4 (IM2I, IM2I, |//|, |//|) 

l;;|2 M2 
+9M/ (IM2I, l/^l) + 41og ^ + Slog ^ + 14 

/*o A'o 

-8Re (MiM*) £'2 (|Mi|, IM2I, |/.|, |/^|) - 8L'4 (|Mi|, IM2I, |/x|, |/x|) 

+iy(|Mi|,|^|) + 3iy(|M2|,|/^|) + 4] 

-12£'2 |Mi|, |Mi|2 - 12D4 (|/x|, |/x|, |Mi|, |Mi|) 



+3iy(|Mi|,|/.|) + 41og 



/^o 



+ 6 



(39) 
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Figure 2: Upper row: A subset of one-loop diagrams for Bq — Bg mixing in the effective two- 
Higgs-doublet model. Lower row: Tree and one-loop diagrams contributing at large tan /3 when 
employing the Lagrangian Cnh and tree-level couplings. The crosses denote the flavor-changing 
neutral Higgs couplings and (in diagrams (f) and (g)) loop-suppressed Higgs mass terms. On the 
lower row, arrows designate the flow of the conserved U{1) charge discussed in Sect. 2.3 



where the loop functions Bq, Cq, Dq, Bq, D2, D2, and W are defined in Appendix B.4, and the 
notation Aj refers to the matching scheme as explained in Sect. 3. Inspecting Eq. (28), A7 enters 
quadratically, which formally is of higher loop order. Nevertheless, it can be seen that \j oc 
as opposed to A2A5 oc g^yf, which can partly offset the additional loop suppression. Indeed we 
find that, numerically, neglecting A7 is not always a good approximation (Sect. 4). 

The form of the matching result depends on the renormalization schemes of both the full 
theory, i.e., the MSSM, and the effective theory, i.e., the 2HDM. The latter cancels in physi- 
cal quantities, while explicit MSSM scheme dependence cancels against the one implicit in the 
MSSM parameters, to ensure that the couplings in the effective theory are independent of the 
renormalization of the MSSM at any given order of perturbation theory. The residual scheme 
dependence in both cases may, however, be important as we are considering a leading effect. We 
will discuss scheme issues in Sect. 3.1, paying special attention to the definitions of tan (3. 



Higgs loops There is a considerable number of one-loop diagrams in the effective 2HDM that 
can contribute to B — B mixing amplitudes (Fig. 2, upper row). These give the following 
contributions to the Wilson coefficients multiplying QY^^ and Q^^^: 
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1 rri n 

^^r^lmggsioops = -4^2eos2 /3(1 +?3tan/3)2 CI^mU^.^^^^^' " ^"^^^ 

In these expressions, we have neglected the small Yukawa coupling and employed tree-level 
MSSM mass relations, in agreement with our approximation of working to leading order in small 
parameters (in the present case, the loop factor l/(167r^)). C^^^ is suppressed by two powers of 
niq/mb inside k*1 in the MFV case, hence beyond our accuracy. The results Eqs. (40) and (41) 
involve a great deal of cancellations, which can be understood in terms of symmetry arguments, 
as explained in Sect. 2.3 below. We note the absence of charged-Higgs contributions in the 
approximation considered here. 



v/Af -suppressed effects All of the couplings given in Eq. (11) correspond to the zeroth order 
in the v/Mqusy expansion, or equivalently to the level of dimension-four operators. Gauge 
invariance forbids dimension-five operators built from quark and Higgs fields, so the leading 
higher-dimensional operators have dimension six. This can lead to more general Higgs-fermion 
couplings than those deriving from Eq. (11) and, in consequence, the cancellation leading to 
Cf^^ — might be broken. To see that this is indeed the case, consider the operator 

Q^'^ - -TT^{HiH^){bRHiQL), (42) 

-'^^SUSY 

which gives rise, inter alia, to effective dimension-three and -four couplings 

1^ bRSL + -^{hnsi^hl + 2 6«s^C)- (43) 

-"^^SUSY -'^-'SUSY 

The first term is removed by a rediagonalization of the quark mass matrices, but the two remain- 
ing terms, in general, are not. The appearance of /i° in addition to hP* leads to a contribution to 
(jSLL proportional to /t^g C'^^\ However, because of i?-parity, SUSY particles do not contribute 
to tree graphs with external standard particles only, such that Q^^^ (or any other higher-dimension 
operator) is only induced at the loop level, and this loop- suppression factor is not compensated 
by factors of tan/?. (Recall that the 0{1) FCNC couplings at dimension four are nothing but 
rotated tree-level Yukawa couplings.) Hence any t'/MsusY corrections that break the cancella- 
tion in JF~ involve an additional loop suppression, and can be neglected for the present analysis. 
On the other hand, as Eq. (43) shows, the higher-dimensional operators do have an impact on 
the rediagonalization of the quark mass matrices and, consequently, on the size of the FCNC 
couplings Kbq. These effects preserve the cancellations in T~ discussed above but have a mild 
impact on the FCNC couplings multiplying in (cf Eq. (35) and the discussion around 
it). 



2.3 ?7(l)pQ and effective Lagrangian for large tan /3 

To better understand the various types of cancellations in and in the Higgs-loop contributions 
to CY^^, as well as the suppression of T^, we now introduce an effective 2HDM Lagrangian at 
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large tan /3. This will allow us, on the basis of simple symmetry arguments, to clarify the role of 
the parameters A5 and Ay, the structure of Eqs. (18), (19), and (28), as well as the vanishing of T~ 

for tree-level Higgs couplings at leading order in 1/ tan (5. It also provides a tool for computing 
loop diagrams involving Higgs bosons efficiently and consistently, which may be useful in other 
contexts such as collider processes with Higgses in the initial or final state. 

by the minimization conditions and trade 



As before, 



we eliminate m^^. 



m 



22' 



and 



12) 



{ml2y for Ml via Eq. (25). We then take the limit 



Vd 



0, v„ — > V, fixed, Xi fixed, 



(44) 



of the Lagrangian (12) in the broken phase."* We also keep the Yukawa couplings fixed when 
considering the couplings to fermions. In this limit we have ^ — H^, ^' — eH^, and 



(45) 



If there were no mixing among neutral Higgses, we would have 0^ = /i° and 0^ = H^, and 
would be a mass eigenstate. The mass matrices are compactly expressed by the quadratic 
potential 



y(2) _ 
^Itb — 



2 . 



+ 



+ \K 
.^2*+h.c. 



A 



2 2 ,2 



(46) 



valid up to corrections of order cos (3 ~ l/tan/3 -C 1. The trilinear terms are given in Appendix 
D; the quartic terms follow trivially from those in the symmetric Lagrangian Eq. (12). Note that 
the first line of Eq. (46) is symmetric under the C/(l) Peccei-Quinn (PQ) transformation 



hi K 



e ^^h^, or equivalently, Ha^e^^Ha, 



(47) 



while the second line is not. In the MSSM, the non-invariant terms appear only at the loop level. 
We note that the C/(l) symmetry is not spontaneously broken in the large-tan /3 limit, so there 
is no massless boson, in agreement with our keeping fixed. ^ Next, a PQ transformation 
makes A5 real, such that the first term on the second line of Eq. (46) contributes with opposite 
sign to the mass terms for 0^ and Xd = —Aq + C(cos splitting the two. There are only two 
independent mixing angles that do not vanish: they can be identified with the CP-conserving 
angle a — C(Ay) and a CP- violating a' — 0{X\); a third angle present in the general 2HDM is 
suppressed by 0(cot /3; v/M). All of these are symmetry-breaking effects. To lowest order in 
the PQ-breaking couplings, the mass matrices are diagonalized by 



( 



V 



M 



Ai.t.2 



'Ml- 



A2t>^ 



1 







1 



<pd 



(48) 



"^This procedure will be justified in Sect. 3.2 

^Also at finite (but large) tan (3, there is no (pseudo-) Goldstone boson, as mf ^ 
mass terms of both and A° (see also Sect. 3.2). 



Ml > contributes to the 
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ml^X2V^, ml = Ml + lX^lv^, ml = M\. (49) 

In a general basis, CP-violating Higgs mixing is present if and only if A7/A5 is complex. Note 
that there is no mixing for the charged scalars according to Eq. (46), i.e. no mixing between 
charged-Higgs and Goldstone bosons due to sparticles in the large-tan (3 limit. 

These considerations can be extended to the Higgs-fermion interactions. The operators up to 
dimension four follow from (11), which, in the limit of infinite tan (5, becomes 

^\th — Mdij HIQlj — dm Qij ■ Ha 

—"^Ri ^Uij QlJ -Hy, + URi KijW^QLj + h.c. . 

This can be made approximately invariant by extending the symmetry transformation (47) to 
fermions. One judicious PQ charge assignment is 

dRi e^^dRi, Qij — > Qlj, ur^ — > ur^, (51) 

which commutes with the SM gauge group, implying that neutral and charged gauge boson 

couplings respect the symmetry. It has been previously used in [13] to classify the Higgs-fermion 
couplings in MFV. However, since for MFV one has one more small parameter Hqb/nhq oc m^/ nib 
for q = s or d, it is useful to consider the following variant of Eq. (51): 

bR e'^bR, qr qr, Qlj Qlj, ur^ URk. (52) 

Now KijdRiQLj ■ Ha in Eq. (50) breaks the symmetry unless dRi — dRs — bR. However, all 
C/(l)pQ breaking is still proportional to one of the small parameters of Eq. (36): Kgj = 0{uj) 
and kij = 0{l). The modified symmetry Eq. (52) forbids all operators in the weak Hamiltonian 
Eq. (4) (Table 1), including the would-be leading one, Qf^^, except for the standard-model oper- 
ator QY^^ and for Qf^^. The last two are, however, forbidden by the original charge assignment 
in Eq. (51). Hence the Wilson coefficients of these operators are suppressed by a; = niq/mb or 
by factors of loop-induced effective couplings, respectively. 

At the tree level (in the 2HDM), JF+, which induces Q2^, is multiplied by a factor k*^, which 
is a PQ-breaking coupling. On the other hand, JF~, which induces Qf^^, is multiplied by the 
unsuppressed factor k^^. Hence must be proportional to PQ-breaking couplings in the Higgs 
potential (up to 1/ tan /3-suppressed terms). This also seen from the fact that in the infinite tan /3 
limit, it is given by 



/ 



d^x(Tf/i°(x)/i^(0))), 



which vanishes if the PQ symmetry is unbroken. Explicitly, in the large tan (3 limit one has: 

j,-^ -(A^A^ - xfy ^ -{X2X1 - xfy 
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Operator [field content] t/(l) charge ^uppression of leading Remark 

^ ^ Higgs-mediated contribution 



Qlf [bRqLbRqL] 


2 


A5 (sparticle loop) 


new 


<5i,2 [bRqLhqR] 


1[0] 




known 


[hqLhqL] 





2HDM loop 


SM operator 


gr^" IbRqRbRqR] 


2[0] 


(jj^x 2HDMloop 


tiny 


QfT hqRbLqR] 


0[-2] 


a;^ X sparticle loop 


tiny 



Table 1: Charges of the operators in the weak Hamiltonian under the approximate U{1) sym- 
metry discussed in the text, see Eq. (52). The number in brackets denotes the charge under the 
"unmodified" charge assignment of Eq. (51). 



where the rightmost expressions hold up to higher orders of small couplings. For JF~, this is 
identical to the sum of the two leading diagrams in a "mass-insertion approximation", where the 
PQ-breaking contributions to the Higgs mass terms are treated as interactions (Fig. 2 (f) and (g)). 

At the loop level (in the 2HDM), up to doubly suppressed contributions one can employ the 
PQ-conserving parts of Eqs. (50) and (46), i.e. set A5 = Ae = Ay = 0, as well as ignore Kgh and 
kij. The matching onto the weak Hamiltonian can be organized according to one-light-particle- 



irreducible chirality amplitudes. There are three amplitudes 

Arr = (T(bR{xi)bR{x2)sL{x3)sLix4)^^, (55) 

Arl = (T(bR{xi)bLix2)sLix3)sR{xi)^^, (56) 

AvLL = (T(bL(xi)bL{x2)sL{x3)sL(x4))), (57) 



plus the parity conjugates of Arr and Avll- (We have omitted amplitudes that cannot match 
onto Lorentz-invariant local dimension-six operators.) Only Avll is invariant under U{1)pq 
(both versions) and can be generated from a symmetric Lagrangian. It matches onto the standard- 
model operator QY^^- There is a single diagram contributing, see Fig. 2 (h). (Diagram (i) 
matches onto Q^^^ and would be allowed for the unmodified PQ assignment of Eq. (51).) 

The present discussion could be extended to other processes and to higher loop orders, by 
systematically treating the PQ-breaking couplings as interactions and working to a fixed total 
order in the small parameters; in practice, at such higher precision, one might want to extend the 
effective 2HDM by higher-dimensional operators to account for v/MSUSY corrections. 

Finally, let us remark that because our choice of shift parameters and minimize the po- 
tential V in the Lagrangian of our effective theory and not necessarily the full effective potential, 
the one-point functions for the (shifted) Higgs fields (0|/ij|0) (hi = (pu, 4>d-, A^) will, in general, 
not vanish. Hence also "tadpole" diagrams involving quark or Higgs loops would have to be 
considered at the outset [Fig. 2 (e)]. That they cancel va B — B mixing in our approximation 
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follows from the fact that no such diagrams are present when working with a complex hP^ field 
and the Lagrangian VJtb- Tadpoles may, however, be relevant in other contexts. We discuss our 
renormalization of Vu, Vd, and tan /3 in detail in the following section. 

3 Systematics of the large-tan (3 MSSM 

The present section is devoted to certain technical aspects of the large tan (3 limit. The first con- 
cerns the definition (i.e. renormalization) of tan in the MSSM and in the effective two-Higgs- 
doublet-model description of low-energy (i.e., Higgs, electroweak, and flavour) phenomenology, 
and the matching between the two. This is of phenomenological importance, as tan j3 defini- 
tions used in the literature on the MSSM are known to differ by parametrically large expressions 
0(tan /3 X loop factor). This can lead to ambiguities in the value of tan /3 of 10-15 in certain 
regions of the MSSM parameter space between schemes that have been extensively used in the 
study of radiative corrections to the MSSM Higgs sector [29]. Having clarified the connection 
between our "full" and "effective" tan/3, we justify the systematic expansion in 1/ tan/3 at the 
Lagrangian level employed in Sect. 2.3. 

3.1 Renormalization of tan (3 

In the MSSM, tan/5 = v^/vd is defined as a ratio of vacuum expectation values. This is an 
unambiguous notion at tree level, because a preferred basis is provided by the chiral Higgs su- 
permultiplets of definite hypercharges ±1/2. Beyond tree level, a scheme dependence arises as 
the bare parameters p° (pj = ■m\,m\,BiJi,g,g',QXc.) are renormalized, p° — Pi+5pi, as well as in 
the normalization of the fields and in defining renormalized shift parameters Vd,Vu. To formalize 
the renormalization program, we first define bare shifts that minimize the bare effective potential 
including radiative corrections, which is equivalent to requiring vanishing one-point functions 
for the shifted fields, i.e., 

{hi - ^v^) = 0, (58) 

such that the are indeed vacuum expectation values. Identifying (for any definition of renor- 
malized shift parameters) 

v^^zl^^ivi-Svi), i^d,u, (59) 

scheme dependence arises through, and only through, field renormalization and the counter- 
terms Svi. Ref. [30] argued that for a stable perturbation expansion it is desirable to define the 
renormalized Vi such as to minimize the renormalized effective potential, i.e. 5vi — 0, and 
implemented this proposal for DR field renormalization and Landau gauge. The same condition 
and gauge fixing was imposed in the computation of one-loop corrections to the MSSM Higgs 
masses in [18-21]. Refs. [22,23] chose to work with on-shell fields and in gauge instead, and 
their shifts do not strictly minimize the one-loop effective potential. In fact, in general gauges, 
for 5vi — the effective action is not finite and the Vi are both divergent [22, 31] and gauge 
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dependent [31, 32] (as are the bare vevs v^).^ Hence to have finite renormalized Vi and tan /3, 
5vi^Q, containing a gauge-dependent divergence, is required. For tan /3, we have 

tan/3° = % tan/3f 1 + Uz^ - - ^ + ^) = tan/3 + aan/3. (60) 

\ 2 2 

Minimal subtraction for Z^, 6vu/vu, Svd/vd defines tan/3°^ [21]. It also follows from 
Eq. (60) that a change between two schemes R and R' can be calculated from 

tan f3^ - tan (3^' = S tan (3^' - S tan (3^, 

hence any scheme where 5 tan (3 is a pure divergence has tan/5 = tan/?°^ regardless of any 
nonminimal field renormalizations as those employed in [29]. In the latter case, however, 6vu, Svd 
are nonminimal and the counterterm for tan /? has no simple relation to the field renormalization 
constants^ 

tan (3^^ is gauge dependent [34], but to one-loop order, the gauge-dependence drops out for 
the R^ gauges. In spite of its gauge dependence, the DR scheme for tan (3 has been shown to 
lead to a well-behaved perturbation expansion [29] and is also used in the most recent version of 
the publicly available computer programs FeynHiggs [35] and CPsuperH [36]. 

A second issue is that a fully minimal subtraction scheme, where in particular = 0, 

generally entails Vi that do not minimize the (renormalized) tree potential, such that the renor- 
malized Lagrangian contains linear terms 

jC D td(t>d + tn(t>u (61) 

for the shifted (real parts of the) Higgs fields. On the other hand, from Eq. (58) and Eq. (59) it 
follows that 

° = + F^ + = (62) 

always holds, if only 5vu and 5vd are included in 5ti. The presence of tu, td is perfectly fine, 
but tadpole diagrams then have to be retained in the calculation. (In particular, they appear in 
the expressions relating Higgs and gauge boson mass parameters to the Lagrangian parameters. 
If all renormalization constants are minimal, Eq. (62) determines U in terms of the bare proper 
one-point functions F*^^^ [21].) Yet it may be more convenient to perform additional finite renor- 
malizations to work in a scheme where ti — 0. This can be achieved either by suitable finite 
terms in 5vi or by finite renormalizations of the mass and coupling parameters. The former shifts 
tan /3 from its DR value according to 

tan P''" = tan /3^ f 1 - ^ + . (63) 

V Vd Vu J 

^This is in particular true in i?^ gauges if ^ 7^ 0. The apparent contradiction to the results in [33], whose authors 
are able to renormalize the effective action with purely "symmetric" counterterms, is resolved by noticing that in 
the Lorenz gauge employed in [33] the gauge-fixed Lagrangian still respects an invariance under constant ("global") 
gauge transformations. This is sufficient to forbid divergences that cannot be removed by symmetric counterterms. 
Conversely, the gauges break also this global invariance, for instance through Goldstone and ghost mass terms, 
which are indeed responsible for the "non- symmetric" divergences at one loop [31]. The exception is the Landau 
gauge ^ = 0, which has the invariance. 
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The latter option does not modify tan /3. 

Going from the MSSM to a general 2HDM, tan (3 becomes - strictly speaking - an ill-defined 
notion as there is no preferred basis. Identifying Hi = —eH^ and H2 = H^, an SU{2) rotation 
Hi — * UijHj removes the vacuum expectation value of one doublet; this corresponds to the 
($, $') basis introduced in Sect. 2. Only $ receives a vev, provides for the Higgs mechanism, 
and has flavour-conserving couplings, while is an ordinary scalar with FCNC couplings. To 
make contact with MSSM phenomenology, however, it is useful to keep the notion of tan /3 in 
the effective theory. In principle, we could fix a basis to enforce tan (3^^'^ = tan P^^, but find 
it technically simpler to allow for a parametrically small (i.e. not tan /^-enhanced) shift, as we 
discuss in the following. 

In complete analogy with the MSSM case dicussed above, if we employ a general gauge and 
MS everywhere in the effective theory, Vu and va will not minimize the tree-level (nor the effec- 
tive) potential. This would require a modification of the formalism in Sect. 2. In particular, in 
writing the mass matrices Eqs. (23-25) and the flavour structure of the scalar- fermion couplings 
in Eq. (1 1) we assumed the minimization conditions h = t2 = 0. To avoid such modifications, as 
well as changed expressions for neutral meson mixing, we can either perform renormalizations 
on the parameters mfi and 777.22 ^^^^ ^^^^ ^1 ^2 minimize the 2HDM potential, or achieve this 
through nonminimal 5vi 2. We pursue the latter option, keeping the symmetric parameters of the 
2HDM minimally subtracted. This has the added virtue that the tan P such defined is gauge in- 
dependent at the order considered, as it is fully determined by MS mass and coupling parameters. 
These are gauge invariant at one loop, which is clear from our explicit matching calculation. We 
presume this to hold also at higher orders, at least if the appropriate wave-function renormaliza- 
tion is employed. The 5vi are determined entirely in terms of "light"-particle loops and, at least 
at one loop, do not lead to parametrically large shifts oc tari^ Pj^, as can be verified from the 
explicit expressions for the tadpoles in [23] or by considering tadpole diagrams in the large-tan P 
effective Lagrangian. 

To find the precise connection between tan P^^ and our effective tan P, consider the total 
tree plus one-loop contribution of the superpartners to the (MSSM) effective action for the gauge 
and Higgs fields. 



^gh 



-i 



;i + ^Zw){-^)W^,W>^''^ + (1 + I^Zb){--)B^.B>^'' 



+{5ij + AZij){D^Hi)''{D^'Hj) - mlHlHj - + ... 

k=l 



■ (64) 



Here Oi are the quartic terms constructed from the Higgs fields appearing in Eq. (12), and the 
dots denote higher-dimensional local terms. The precise values for the coefficients depend on 
the MSSM renormalization scheme. We assume the MSSM has been regularized by dimensional 
reduction while the Higgs fields and tan P are minimally subtracted (DR). The corresponding 
expressions are reported in Eq. (21) (tree level) and in Appendices B.2 and B.3 (one loop). 

Eq. (64) can be identified with the classical action (ignoring 2HDM loops) for an effec- 
tive two-Higgs-doublet model with noncanonically normalized fields. To obtain from this the 
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MS-renormalized Lagrangian in the presence of light-particle loops one simply has to add the 
contributions (which are local) due to loops of 2e scalars present in DRED ^ and subsequently 
rescale the fields. 




(65) 



subject to the condition Z^(l + AZ)Z — 1. This provides the relation between the DR fields of 
the MSSM and one out of an infinite choice of MS fields in the effective theory, labeled "eff". 
We fix the freedom to choose the Higgs basis in the effective theory by setting Zdu = and 
= ^dd ^ 0- The relation between the shifts and tan (3 of the MSSM and of the 2HDM are 
now determined according to 
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v^{^) = vM'^ = - SZddVd + Svl^ 

tanPiiiY'' = tan/^D^ 1 ^— + ^- + SZdd - 5Z^^ - 5Z^d cot {3 



(66) 



Here we have expanded Zuu/dd = 1 + 5Zuu/dd and Zud = ^Zud, and the 6Zij are related to 
the AZij via AZii/22 = —'^^Zuu/dd and AZ12 — —5Z^, with the explicit expressions given in 
Appendix B.l. The shifts Svl^ are defined implicitly as discussed above. In summary, we have 
constructed a tan P which is appropriate for effective weak interactions, gauge-independent and, 
up to an ordinary (i.e., not tan /9-enhanced) loop correction, coincides with the widely used 
tan It means that the tan /? measured in flavour physics, for instance through BR(B^ — > 
IJi^lJi~), and employed in our analysis, can be identified with the corresponding DCPR parameter 
at large tan /3, up to small corrections. 

We note that our framework leads to a transparent expression for the relation between the DR 
scheme and the so-called DCPR scheme employed in [22,23] in the limit v <^ Msusy- In the lat- 
ter scheme, finite but, unlike in our effective 2HDM, "diagonal" wave function renormalizations 
of Hu, Hd are performed, i.e., in our notation, SZ^d — SZdu — 0. Moreover, the renormalization 
conditions include 

— = — , ReEAozo(Ml) = 0, (67) 

Vu Vd 

where 2 40^0 (/c^) parameterizes the mixing according to E^02o(^) = A;'* E^02o(/c^). Now, 

the sparticle contribution to S^02o(A;^) reads 

Saozo (A;') = sin' /3 AZdu + sin /3 cos - SZdd) + ■ ■ ■ , (68) 

where the dots denote terms proportional to cos (3 but not involving the wave-function renor- 
malization constants. This follows either by considering the mixed gauge boson-Higgs boson 



^Integrating over the 2e scalars leaves a path integral over light fields that is identical to that in the DREG- 
regularized effective theory, including the 1/e divergence structure. We recall that the 2e scalars should be thought 
of as having a nonzero mass of 0(Msusy) [37]. 
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bilinear terms resulting from the covariant kinetic operator for the Higgs fields in Eq. (64), or via 
the Ward identity 

A;^E^o^o(A;') + Mz^a^Aoik^) = 0{e - Ml) (69) 

(which is trivially satisifed in our SU (2)-invariant formalism) from the terms bilinear in the 
gauge fields in the same term. The two conditions in Eq. (67) then determine 



5 tan = tan - 5Zdd) = - tan^ (3 Re AZd„ + . . . , (70) 

where the omitted terms are not tan /9-enhanced. This explains the large numerical differences 
between tan f-P^ and tan fp^'^'^ found in [29] as a parametrically large effect. Hence, tan (3 
measured in flavour physics should not be identified with the corresponding DCPR parameter at 
large tan /3. 

As with the Higgs fields, we explicitly decouple the contributions of heavy particles to the 
gauge field wave functions (hence to g^'^) by a finite renormalization g^'^^ — g^'^ + g^'^5g^'\ can- 
celling the terms AZb and AZw in Eq. (64) of the gauge fields, = and = Z^W^. 
For DR-subtracted MSSM couplings, this gives MS-renormalized 2HDM gauge couplings. 

We denote the quartic couplings in our 2HDM scheme by Aj. The finite renormalizations 
leave A5 invariant, A5 = A5, while the other quartic coupling constants transform like 

Ai = Ai + g'SZ^,, + 1 (^g% + g'^Sg') , A2 = A2 + fSZ:, + \ (gHg + g'^Sg^ , 

A3 = A3 - ^ {SZ^^ + SZ:J - \ {g^g + g''5g') , A4 = A4 + / + + gHg, 

Ae = Ae - -^^Z*^, A7 = A7 + —bZ^d , (71) 

where x'^ and denote the real and imaginary part of x respectively. The couplings Aj are MS 
couplings from the viewpoint of the effective theory. 

The modification of the dimensionless couplings by the finite wave function renormalizations 
affects the B — B mixing amplitudes as a formally higher-order effect, as does the scheme 
dependence of tan 13. Unlike the latter, however, the former is never tan j3 enhanced unless the 
wave function renormalization constants themselves are. 

Invariance of B — B mixing under field renormalization 

The effects of Eq. (65) on the Higgs-mediated FCNC Eq. (1 1) are twofold: (i) the values for cos (3 
and sin /3 in Eq. (10) are modified. This cancels the contributions to from the redefinition of 
the mass matrices up to a global factor: 



J^+(A,, v^,d. Ma) ^ ^+(A:, <<i, M'j,) = — |det Z\' v^,d, Ma), (72) 

.F-(A„^;„,d,M^) ^^-(A:,<„M;) = ^{detZ*)'T-{X,,Vu,d,MA). (73) 
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(ii) the Vu,d renormalization in (i) comes with a modification of Kij: 

t^ij^<j^'^idetZ*)-'i^ij- (74) 

The above factors cancel each other out in the products T~ and k*^ Kf,, , as they should. 
In particular, our choice of wave-function renormalization acting on the leading FCNC coupling 
Eq. (29) produces an extra term: 

8K^q = -Kfeg (sj - SfiCfi 5Zl^ + 4 ^Ku) ■ (75) 

Considering Eqs. (29), (30), (71), and (75) gives the same Wilson coefficients and Cf^^ 
as does considering Eqs. (29), (30), and (71) with the finite parts of 5Zij set to zero. While in 
practice, wave-function renormalization has to be performed to relate the parameter Ma to the 

— 1/2 

physical Higgs boson masses and to take v — {\/2Gf) — 246 GeV beyond leading-order 
precision, such renormalizations are not the source of a non- vanishing of the amplitude, to 

be found instead in the corrections to Higgs masses and mixings (via the self-couplings Aj, in 
particular A5); wave-function-renormalization effects enter that amplitude only at higher orders 
(as might have been expected). In this regard our findings disagree with the conclusions of [25]. 



3.2 Health of the large-tan (5 limit and fine-tuning 

In Sect. 2.3 we took the limit tan (3 ^ 00 (v^ 0, M\ = const, + v"^ = const, Aj = const, 
Vu and Vd defined as minima of the tree potential) at the Lagrangian level. One might wonder 
whether this procedure is valid at the quantum level. To justify it, we show that the Vd — case 
and the Vd case are analytically connected, i.e. one can be reached from the other without 
a phase transition. It then follows that amplitudes are (in some neighbourhood of a parameter 
point with = 0) analytic functions of the parameters (either "symmetric" or "broken"). The 
renormalizability of the effective Lagrangian Vitb then follows by standard arguments from the 
fact that it is equivalent to the symmetric Lagrangian Eq. (12) (for a certain choice of parameters), 
which is renormaUzable. 

We first note that the number of independent minimization conditions is unchanged in the 
Vd = limit. First, for general values of the parameters, out of the four real (two complex) 
minimization conditions, at most three are independent. This follows from the U{1)y invariance 
but is easy to verify explicitly. Fixing Vu to be real and positive, three polynomials of degree 
three determining three unknowns Vu, v'^, remain. The system has a solution Vd — Oif 

X,ml, + Xrml,^0, vl^-^-^. (76) 

Here the second equation determines w„ = w as a function of and A2 similarly to the case of 
a single doublet, while the first equation can be viewed as a fine-tuning condition between 
and Ay. The dimensionless, complex parameter 

m?2 A7 mno 

e = ^ + T^— §^ (77) 
mfi A2 mfi 
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parameterizes the deviation from the fine-tuning Umit; we may trade in favour of e. Clearly, 
at the limiting point e = we indeed have three independent equations. Now, it is easy to verify 
that, writing the four real minimization conditions in the form 

fi{\i, m^ij, Vu, Vd) = 0, (78) 

the Jacobian matrix 



d(fl,f2,f3,U) 



(79) 



has maximal rank (3) at any point with Vd — 0. (Physically, this just means that the neutral Higgs 
mass matrix has three nonzero eigenvalues there.) Hence, by the implicit function theorem, 
we may solve for (?'„, Vd) in a neighbourhood of it, where the solutions will be (real-)analytic 
functions of e. In particular, Vu behaves analytically (and is strictly positive) around e = 0, i.e. 
no phase boundary is encountered. Explicitly and to linear order, the real and imaginary parts of 
Vd are determined by 



1 + ^ 



^ + 2 m^i 




+ 0{e'), (80) 



such that tan P — 0{l/\e\). The nonsingular linear term allows us to change variables from 771^2 
to a complex Vd. Of course, we may always perform a field redefinition of Hd such that Vd is real. 
Then, the mass parameters besides ml^ are power series in 1/ tan P, which read 

2 _ VdUlli I 1„, /\ I \r\ I ^„.2\r , /nr..2 /„.2 



^22 = 



+ ::VuVd (A3 + A^) + -viX^r + 0{v'dK), 



2 ^ 0/ ■ 2 

2 _2 , A3 - A5 2 



Mi = rni, + + 0{vdlvu), ^^^^ 



Ml,^m\, + ^^^±^v' + 0{vd/v:). 

We see explicitly that we can continuously change the dimensionful parameters in the Higgs 
potential from a situation where Vd^^to one where = 0, keeping M\ (and the dimensionless 
couplings) fixed, as was assumed in Sect. 2.3. The last three equations illustrate that the large- 
tan (3 case is characterized by a "primary" doublet which receives a large vev Vu and a 
"secondary" doublet Hd with a positive gauge-invariant mass that receives corrections of 
0{v'^) and 0{e), respectively, due to its dimensionless and dimensionful couplings to H^. Those 
corrections differ among the physical components of the doublet, approximately to be identified 
with H^, A^, H^, due to electroweak symmetry breaking. In principle, mf^ could be negative, 
but in that case, t'd ~ will typically not be the global minimum of the potential. 
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We close this section by considering the fine-tuning which is necessary to obtain a large tan (3 
while keeping the mass Ma fixed. For real, Eq. (80) implies 

m{2 = -^^"^22 + cot/3(mii H vj, (82) 

which illustrates the tuning that is known to be necessary to have large tan P in the MSSM. For 
the generic situation mf ^ ~ ^susy ^ the right-hand side is dominated by the m^^ term: 
A7 is down by a loop factor relative to A2, and 777,22 ^ v'^ ^ -^susy (the little hierarchy). Hence, 

^i27^sW--l/tan/3, (83) 

which implies an extra tuning beyond the one to achieve the correct weak scale. For smaller 
777^j^ ~ M\ ~ M|, which is interesting from the point of view of i?-physics phenomenology, the 
required tuning gets even worse - unless, of course, the whole SUSY scale is lowered to the weak 
scale, which is, however, problematic since then X2v'^ is generally below the experimental 

lower limit. On the other hand, as we have seen, ~ 777^^, such that no extra tuning is required 
to keep Mj^ finite, while one might have expected otherwise from the well-known tree-level 
formula 

Ma = (tan P + cot P)mj2 , (84) 

which is generalized by Eq. (25). Also, while a small 777^2 indeed sensitive to radiative correc- 
tions, those are automatically correlated with shifts of and in consequence of tan (3 in such a 
way that M^ receives only mild corrections. 



4 Phenomenology 

In Sect. 2, we performed a detailed study of the supersymmetric contributions to AM^ and 

A Ms in the generic framework of an effective 2HDM. The corresponding matching coefficients 
were computed at the one-loop level in Sect. 3 and Appendix B. In this section, we assess the 
maximal size of the various types of effects identified in the MFV case taking into account the 
existing constraints on the supersymmetric parameter space, in particular from the Bg — > H'^jJT, 
— > t^Vt, and h ^ branching fractions. For convenience, we start with a compendium of 
the formulas derived in Sect. 2.2: 

AMg = I AMgS^ + AMj^^ + AMj^^ + AMf ^ | = 1 1 + /i^ | l^Mf^ , (85) 
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where the standard-model, the left-right and left-left Higgs-pole, and the neutral Higgs-loop 
contributions read 



VLL 
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246 GeV. The flavour-changing and flavour- 



conserving quark-Higgs couplings were defined in Sects. 2. 1 and 2.2: 
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VCp 1 + 63^/3 ' 
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with eo,y and ?3 given in Eqs. (3 1-33) and yt in Appendix A. The J^^ factors describing the prop- 
agation of the neutral Higgses were defined in Eqs. (18) and (19), with the effective couplings Aj 
entering the neutral Higgs mass matrix computed in Sect. 3.1 and Appendix B. For large tan /3, 
we have in very good approximation: 



Ml' 



Mi 



(88) 



(exact formulas were used in our numerical analysis though). Explicit expressions for A5, A7, 
and A2 in the MFV case were given in Eqs. (39-37). Altogether, counting ey , T~ ~ (IGtt^)"^ 
and Ma ~ 120 GeV to get an idea of the naive size of the various effects in the absence of 
constraints, we obtain: 



2.40 



0.053/2.75 



VLL 



3.2/0.71 



+ 0.35 



167r2(-A^ + A*7VA2)(120 GeV)^e^''l^ 



Ml 



pSLL /pVLL 

-1.36/0.71 
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(1+63%)^ 
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(89) 



rub 



1 2 



2.75 



where is in GeV and ^ = arg(7«). hd is given by the same expression with replaced by 
rrid, so that the first term becomes subleading. 

A first obvious remark is that AM^'" cannot compete with AM^^^ or AM^'" unless yb be- 
comes non-perturbative. This is rather accidental (notice the small loop factor in Eq. (86) as well 
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as the smallness of P^^^ with respect to and ^^). Further, the contribution of AM^^ 
seems somewhat limited. However, the loop functions A5 and A7 could be enhanced for large // 

or 5, see Eqs. (38,37). A more quantitative analysis is thus desirable. In the next two sections, 
we perform a random scan of the MFV-MSSM parameter space to find the maximal AM^^ and 
AM^^ values allowed by current experimental data. Eqs. (85-89) do allow for new CP-violating 
phases^, yet these will be set to zero in the scan. CP- violating effects within the MFV scenario 
will be shortly discussed in Sect. 4.3. 

4.1 Scan of the parameter space 

The values of the various input parameters used in the scan are collected in Table 2. Note that 
only the products Pj"^^* ^^'^ P2^'rnd, or alternatively P2^ms and md/nis, are needed, see 
Eq. (86). We scan over P2^ms but keep nid/ms fixed as AMj^ is doomed to be small any- 
way. The decay constants and CKM factors iVt^T^^I are not specified. Instead, outputs are 
formulated in terms of ratios free from these rather poorly known parameters. Finally, we take 
a = 1/127.9, sin^ Ow = 0.231, and Mz = 91.1876 GeV. 

For simplicity, the gaugino mass parameters are assumed to have the same sign (which we 
can choose positive), as well as the trilinear terms (positive or negative). Note that the absolute 
scale of MsusY plays no role as supersymmetric parameters enter eg y and A; by means of ratios. 
Only the spread of the interval chosen for Msusy matters. Still, Msusy should not be taken too 
large to help satisfy the b 57 constraint in the case /x < 0. We will come back to this point 
later. We allow for rather large values of Ma, close to the lower end of the interval chosen for 
^susY- Still, the matching performed in Sect. 3 and Appendix B remains valid as the corrections 
from higher dimension operators at the electroweak scale are ruled by the ratio w /Msusy and not 
Ma/ Msusy- The formulas for the various observables at the B mass scale are thus unaffected. 



Quark masses and as 


Bag factors 


SUSY parameters 


mt = 164 GeV 
rub = 2.75 GeV 
rrid/ms = 1/19 
as = 0.108 


P^'^^ e [0.66,0.76] 
P^^nis e [0.12,0.22] GeV 
pSLL g [_i.48,-l.24] 


tan/5 e [10,60] 

Ma e [120, 600] GeV 

Msusy e [600, 1800] GeV 



Table 2: Input values. Here Msusy stands for any of the supersymmetric parameters M^^, 
Mj^, M-j^^, M^j^, M^j^, \at\, \ai)\, |a^|. Mi, M2, M3. As explained in Sect. 3.1, the renormalized 
parameters Ma and tan /3 are identical in the MSSM and effective 2HDM for Msusy > v. The 
quark masses and are defined in the MS scheme at the scale m^. The bag factors, defined at 
the scale mt as well, are discussed in Appendix C. 

The constraints imposed on the points generated inside the above ranges are summarized in 
Table 3. We now discuss them in turn: 



^Let us recall that in that case M\, defined as the non-zero eigenvalue of the CP-odd mass matrix Mf in Eq.(14), 
is no more an eigenvalue of the full Higgs mass matrix. 
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i) The bottom Yukawa coupling yt, is maintained small enough, say, yi, < 2, to guarantee the 
validity of perturbation theory. This condition removes possible fine-tuned points in parameter 
space for which the denominators in Eq. (87) are close to zero. 

ii) The lightest Higgs boson mass Mh has to come out large enough to comply with the LEP 
II experimental lower bound. is obtained from the CP-even Higgs mass matrix in Eq. (23), 
with the effective couplings Aj computed at the one-loop level. Higher order corrections to 
A2 are known to be important [18-20]. However, h° comes up in the FCNC vertices Hij of 
Eq. (11) along with a cot /3 suppression factor. The tan /3-enhanced effects considered here are 
thus largely uncorrelated with Mh. For this reason we do not correct the one-loop formulas and 
simply impose > 115 GeV. 

iii) The following bounds are imposed on at and to avoid the occurence of color symmetry- 
breaking vacua at tree-level [39]: 



(90) 



\at\^ <3iMl + Ml +m:^^ 
lobl^ < 3(M|_ + + m^i 

The corresponding bound for is not imposed as sleptonic parameters have very little impact 
on the quark FCNC considered here anyway. 

iv) The most stringent constraint on the FCNC coupling k comes from the Bs — > A*"*"//" 
branching fraction, which we normalize to AM^ to avoid the occurence of the parameters fs^ 
and VtsVfl. This time the Higgs-pole contribution overcomes the standard-model and Higgs-loop 
pieces. In addition these last two interfere destructively, so we will neglect them. The counterpart 
of Eq. (86) then reads (with m^/M|^ = for simplicity): 



= R,B{B,^l,+l,-f^ 
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where J^p and J^s refer to the Wilson coefficients of the effective operators Qp — {bRSL){i^5i) 
and Qs — (J)RSL){i() arising from neutral Higgs exchanges and 
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167r2 Ml 



2 \|Mi|2'|Mi|2 




H2 -TT7f^i^^2 . (92) 



IGtt^Mi I |Mi|2' |Mi|2 / 327r2M2 llMsP'lMsP 



This result agrees with [40] but disagrees with [41]. The loop function H2 was defined in Eq.(34) 
and Mj2i(j^) = M^iiR) in our MFV scenario. In the large tan P limit and at tree-level in the 
Higgs potential, we have: jFp = — JF5. = J^~^/2 = 1/M\, so that B{Bq /U+yu") is tightly 
correlated with AM^^ [5]. Going beyond the tree-level and large tan (3 approximations we 
obtain: Tp — sp{T'^ — !F~)/2, with given in Eqs. (18) and (19). This formula is actually 
valid in any 2HDM, including arbitrary CP-violating phases (in the CP-conserving case it reduces 
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to the usual identity Tp = Sj^/Mj^. We did not find such a general and simple form for yet 
it is straightforward to write it in terms of Ma, tan/3, and the Aj's (alternatively one can of 

course express it in terms of the neutral Higgs masses and mixing angles). Note that one still has 
J^s = —[T^^T~)I2 up to cot /j-suppressed terms. Sparticle loop corrections to the Higgs self- 
energies turn out to be relevant in the case of Tg- they can be as large as 15% for small Ma after 
all constraints are taken into account, as we will see in Sect. 4.2. Numerically, Higgs-mediated 
effects can easily be very large: 



R, = 9930 



\ ^ (-^5 + IAtIVAs)^^] |167r2ey|' (120 GeV)^ 



Ml 



l+?3t/3|'|l + eot/3|'|l + eM^/3|'Ml 
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40 



(93) 



and ~ Rg. The first correction factor above captures the bulk of the effects from the Higgs 
self-energies, yet the exact formula for Ts should be used for better precision. In practice, the 
looser constraint E{Bs ^'^ ^^)//\Ms < 5.7 x 10"^ ps, obtained from B{Bs fi'^fi'^'^ < 
10~^ [42] and AMfP = 17.77 ± 0.1 ± 0.07^5"^ [43], is built-in in the scan procedure, then 
the current 95% C.L. bound B{Bs ii+ii-)/AMs < 3.3 x lO"'^ ps corresponding to B(Bs 
H+li^y'P < 5.8 X 10-^ [44] is imposed. We also checked the bound B{Bd ii+ii~)/AMd < 
3.6 X 10-^ ps, corresponding to B{Bd ji+ii-y''^ < 1.8 x 10"^ [44] and AM^^ = 0.507 ± 
0.005ps^^ [45]. This provides no additional constraint. Neither do B{Bs^d ~^ 1^^ ^iid AMg^d 
taken separately due to the large parametric uncertainties from /^^ . 

v) The b ^ s'y branching fraction with the energy cut > 1.6 GeV is computed using 
the fortran code SusyBSG [46]. Higgs-mediated effects now appear at loop-level with smaller 
powers of tan /3, so that purely supersymmetric loop corrections (scaling as 1/Msusy) are com- 
paratively more important. For at/i < and relatively light Msusy^ chargino and charged- 
Higgs loops can interfere destructively and more room is left for New Physics. This interplay 
is welcome when ^ < as the charged Higgs contribution then tends to overshoot the ex- 
perimental branching fraction. On the other hand, in that case, the discrepancy between the 
(g — 2)^ standard-model prediction and its present measurement [47] increases (for a recent 
discussion, see e.g. [40] and references therein). The significance of this discrepancy, how- 
ever, is questioned by the new e+e^ 7r+7r~7 BABAR data [48]. We therefore still include 
the situation < in our considerations. The B{b — > 57) experimental world average reads: 
B{b 57)^"? = (3.52 ± 0.23 ± 0.09) x 10"^ [49]. The standard-model central value of the 
SusyBSG program agrees well with the next-to-next-to-leading order prediction B{b — > 57) = 
(3.15 =b 0.23) X 10""^ [50]. We combine the experimental error with the uncertainties discussed in 
Ref. [46] and obtain the following two-sigma range: 2.71 x 10"^ < B{b S7) < 4.33 x 10""^. 

vi) The B~^ — > t^Ut branching fraction is given by 

B(B+ ^ r+i..) = ^TB^ fl.MB+ml {l - |1 - gp\\ (94) 



where 



9P = , . \ ^..2 (95) 



(1 + e^^tp) (1 + ertp) M\ 



2 

H+ 



31 



parametrizes Higgs-mediated effects, e,- is obtained from in Eq. (92) by the replacement 
MfiL{R) — ^ MtL{r)- Corrections to the Higgs potential merely change the value of Mh+, which 
becomes a function of Ma, tan (3, and the various supersymmetric parameters. Again, we in- 
clude these corrections in our numerical analysis. Given the large theoretical and experimental 
uncertainties, we impose: gp < 0.36 U 1.64 < gp < 2.73. The constraint from B{B — > Drv) 
allows to reduce the second interval, and we end up with gp < 0.36 U 1.64 < gp < 1.79 [51]. 



Built-in constraints 


Additional constraints 


B{Bs fi+fi-)/AM, < 5.7 X 10->s 

yfe < 2 

M;, > 115 GeV 

Stability bounds, see Eq. (90) 


B{Bs ii+fi-)/AMs < 3.3 X 10-^ ps 
2.71 X 10-4 < B{b S7) < 4.33 x lO'^ 
gp < 0.36 or 1.64 < gp < 1.79 



Table 3: Constraints built-in in the scan procedure (left) and imposed afterwards (right). 



4.2 Size of the new contributions 

The various Higgs-mediated contributions AM^^, AM^^ and AM^^ normalized to the Standard- 
Model prediction AM^^ are displayed in Fig. 3(a-c) as a function of the FCNC coupling k. As 
expected from Eq. (89), Higgs-loop effects are very small (the bottom Yukawa coupling actually 
does not reach its upper bound yj, — 2in the presence of the other constraints, see Fig. 3(d). The 
upper and lower branches correspond to < and /i > 0, respectively). Further, the contribution 
of AM^^ appears to be much smaller than that of AM^^ despite the fact that mi,T~ can compete 
with rrisT^ , see Fig. 3(e,f). This suppression is a consequence of the Bg — > iT constraint. 
Indeed, large values of T~ are obtained for small values of Ma, to which BiB^ l^'^ l^^) is 
particularly sensitive. As a result, the recent CDF bound [44] only leaves room for very small 
7t couplings, killing practically all effects in AM]^ (and actually also in AAf^^ for such small 
Ma values). In Fig. 3(g), we illustrate this decrease of the maximal k value allowed by the 
Bg — > A*^A*~ constraint with Ma- Blue/magenta/red (dark grey / light grey / grey) points cor- 
respond to Ma — 550/350/150 GeV (the constraints in the right column of Table 3 were not 
imposed to keep the focus on Bg — > A*"*"/^"). As one can see, for Ma fixed, the largest possible 
K first increases with tan/?^, as expected from Eq. (87), saturates the B{Bs — > /^"""/x") exper- 
imental upper bound for some tan /5 value, and is then forced to decrease. For smaller Ma, 
the Bg — > constraint is more stringent and only a smaller l^max can be achieved. This 

growth of 'Kmax with Ma is sufficient to overcome the 1/M| suppression factor in AM\^ but 
not the 1/M| one in AM^. Overall, Higgs-mediated effects in AM^ are of the LR type and the 
room for such effects increases with Ma- The correlation between AAf^ and B{Bg — > /x"*"/^") 
pointed out in Refs. [5] is thus preserved, up to the relatively small Higgs self-energy corrections 
to Bg — > \Ji^ ijT mentioned above Eq. (93). These are only relevant for > 0, tan/5 < 25, 
small Ma, and large A5, though (see Fig. 4). The mass difference in the B^, system, on the 
other hand, remains unaffected. These results seem to contradict those of Ref. [52], where large 
LL-type effects were claimed. Without attempting a close numerical comparison (the sign of 
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AMl"^ / AMl'^ in [52] is actually reversed with respect to ours), let us point out that, as shown 
by Figs. 3(e,f), a large AM^^'^/AM^^ ratio does not automatically lead to large non-standard 
effects in AMg due to the B{Bs — > /^"""yu") constraint. 

Being of the LR type, the maximal effect allowed in AMg is essentially determined by the 
current B{Bs l^'^l^') experimental upper bound for a fixed (but large enough) value of the 
ratio tan/3/M^. This is illustrated in Fig. 3(h) for a slightly larger bound (cf. left-hand side 
column in Table 3). The correlation itself is displayed in Fig. 5, where each diagonal corresponds 
to a fixed value of the ratio ian (3 /Ma- We distinguish the cases > and /i < as the 
latter leads to larger effects due to smaller denominators in Eq. (87) but is disfavoured by the 
measurement of [g — 2)^, as mentioned previously. The sign of the various a-terms, on the 
other hand, has only little impact. Still, in the case < 0, > helps satisfy the 6 — > 57 
constraint. Note that the effect of the — > t^v,,- constraint is particularly transparent on Fig. 5: 
it removes the points with large tan/3/M4 ratios, i.e., the steepest diagonals. Altogether, for 
Ma < 600 GeV, Higgs-mediated effects in AM, can reach ~ 7% (~ 20%) for > (^u < 0). 
These findings agree with those of Ref. [53]. They merely follow from the B{Ba — > iJi^n') 
constraint, as one can see from Figs. 3(g,h) and 5. 

Finally, for completeness (or out of curiosity), we display in Fig. 6 the dependence of various 
quantities on effective couplings or supersymmetric parameters. In particular, in the last four 
plots, we illustrate how the loop functions ^y^ ^t, and A5 increase with the range chosen for 
-^susY (more precisely, they increase with the trilinear and /x terms and decrease with the squark 
and slepton mass parameters Mj^ and Mj^ with f = t,b, r). 

4.3 CP-violating effects 

The Higgs-mediated B—B mixing amplitudes studied here can in principle generate new contri- 
butions to the CP- violating phases measured in the Ba J/ipKs time-dependent CP asymme- 
try and the Bs J/i}(t> time-dependent angular distribution. The coefficients of the sin(AMg t) 
terms are 

Sj/^Ks = sin(2/3 + ^^^^ 
Sj/i,<f> = - sm{-2(3s + 0f ), 

where P = arg[-(y,*y,,)/(K*K,)], /3, = - 8.rg[-{V*Vt,)/{V:,V,,)], and 

0f = argiMfjMfi"^) = arg(l + %). (97) 

In Bg — > J/ an angular analysis separates the different CP components, the sign quoted for 
Sj/ip4, in Eq. (96) refers to the dominant CP-even component. These phases have received a 
lot of attention recently. In particular, the new measurements of —2/3, + 0f by the CDF and 
DO collaborations [54], both more than 1.5 sigma above its SM prediction [55], have triggered 
speculations about the validity of the SM [56]. A possible tension between the value of sin 2/3 
obtained from Sj/^Ks and the amount of CP violation in the kaon system was also pointed 
out [57]. 
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Figure 3: Study of Higgs-mediated contributions to AMg (see text). Black dots denote the points 
in parameter space that satisfy all constraints, while grey dots refer to those that only satisfy 
the initial constraints (see Table 3). In plots (g) and (h), blue/magenta/red (dark grey / light 
grey / grey) points correspond to Ma = 550/350/150 GeV, respectively. Plain lines indicate 
the B{Bs /^^A* ) constraint. In plot (h), the dashed line corresponds to the more stringent 
B{Bs — > /i^yU^) constraint in the right column of Table 3. 
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Figure 4: Correction factor to B{Bg fJ'^fJ'^) arising from supersymmetric loop effects in 
the Higgs self-energies as a function of the effective coupling A5 (left) and the total Higgs- 

mediated effects in Bg /i+yU" (right), for fx > 0. On the left-hand side, the upper (lower) line 
corresponds to Ma = 120 (600) GeV, assuming the approximate formula of Eq. (93) for the 
correction factor. 
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Figure 5: Correlation between AMg and B{Bq — > ): (a) q = s, 11 < 0; (h) q = s, n > 0; 
(c) q — d, fjL < 0; (d) q — d, ijl > 0. The descending lines correspond to a fixed value of 
theratio tan/5/MA. From left to right: tan/^/M^ = 0.05,0.075,0.10,0.13,0.21 GeV-\ The 
ascending lines refer to the B{Bs / AM^ constraints, see Table 3. These lines do not 

take into account the uncertainties on the quark masses and bag factors, nor the effects from 
sparticle loop corrections to the Higgs potential in AMg, B{Bs — > fx^n'), and to the lepton 
Yukawa couplings in B{Bs — > //"*"//"), so that the actual points do not follow them exactly but 
are somewhat scattered. 
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Figure 6: Dependence of various quantities on effective couplings or supersymmetric parameters. 
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Looking back at Eqs. (17-19), it is clear that the new phases cf)^, when associated with the 
Q^^ effective operator, have to be brought up by the quark-Higgs couplings Kij as ^+ cannot 
develop an imaginary part. When associated with Qf^^ or Qi^^, on the other hand, they can 
arise from both the Yukawa sector and the Higgs potential via Within MFV, K*i,K,bq = 
\R\'^\'^fjmqmb/v'^ and only Cf^^ can produce a new phase (via eo,y or X^j). However, the Bg — > 
fx^fj,~ branching fraction is barely affected by CP-violating effects, so that its constraints on 
\k,\ are still very well approximated by the plain lines in Fig. 3(g) (for some representative Ma 
values). As a result, just like in the CP-conserving case, the net effect of the suppression of \R,\ 
and enhancement of J^^ for small Ma is quite small. The MSSM with large tan (3 and MFV is 
thus not able to account for a large non-standard phase in Bg — Bg (or Bd — Bd) mixing, if the 
evidence for such a phase were confirmed. Let us emphasize, however, that the formulation of 
MFV adopted here does not coincide exactly with the full symmetry-based definition of Ref. [13]. 
In the formalism of Ref. [13], it was shown recently that new phases could appear in the 5^^^^ 
sector, in addition to those in the (5^/?)^^'^^ sector [16]. The possible impact of these MFV 
phases via k*;, Kbq in ^ is a priori rather limited due to the Bg A*"*"/^" constraint, yet a more 
quantitative analysis is desirable. 

Beyond MFV, the (^2^ contribution is expected to dominate. As said before, supersymmetric 
loop corrections to the Higgs propagator JF+ do not bring in any new phases. These can only enter 
via the quark-Higgs couplings /t*^ and Kbg. The possible size of CP- violating effects generated 
in this way without violating the existing constraints deserves a study on its own. We will not 
discuss this further here. 



5 Conclusions 

We have studied supersymmetric loop corrections to the MSSM Higgs sector. While the tree- 
level Higgs sector of the MSSM is a 2HDM of type II, the soft supersymmetry-breaking terms 
lead to new loop-induced couplings which result in a generic 2HDM with FCNC couplings of 
neutral Higgs bosons to quarks, even if the supersymmetry-breaking sector is minimally flavour- 
violating. The strength of these couplings grows with tan /3 and precision observables of flavour 
physics are known to severely constrain large-tan (3 scenarios of the MSSM. The appropriate 
tool for such studies is an effective Lagrangian, which is derived by integrating out the heavy 
supersymmetric particles. The abundant literature on the subject has primarily focused on the 
flavour-changing Yukawa couplings [1-7]. Among the FCNC quantities, B — B mixing plays a 
special role, because the apparently dominant contribution of Fig. 1 vanishes. Therefore B — B 
mixing is sensitive to subleading effects, whose systematic study was the main motivation for 
this paper. Pursuing this goal we have derived several conceptual and analytic results which can 
be applied well beyond this topic. They can be classified into three categories: 

1. MSSM Higgs sector 

We have matched the complete MSSM Higgs sector, i.e. both the Yukawa interactions and the 
Higgs potential, onto an effective 2HDM. Our results for the effective Yukawa couplings are 
valid for arbitrary CP phases of at, and the gaugino masses; and Eqs. (31) and (32) correct the 
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gaugino contributions to eo and ey quoted in Ref. [25]. The complete one-loop matching correc- 
tions for the quartic Higgs couplings for the most general MSSM are explicitly listed in one place 

for the first time. This result goes beyond minimal flavour violation and beyond the large-tan (3 
limit. It is well-known that improper choices of the MSSM renormalization scheme can lead to 
radiative corrections which grow with tan /? rendering perturbative results unreliable [29]. At the 
heart of this problem is the feature that tan /3 is an ill-defined parameter in the general 2HDM, 
which permits arbitrary rotations among the two Higgs doublets. In the matching of the MSSM 
onto the 2HDM this feature enters through the wave function renormalization, and we propose an 
MS renormalization of tan f] in the 2HDM which is stable in the limit of large tan (3. The relation 
to a DR-renormalized tan (3 in the MSSM is discussed including electroweak corrections. We 
identify the places in the effective Higgs potential where physical tan /3-enhanced effects occur. 
The coefficients A2, A5 and Ay, which are important for B — B mixing, are explicitly specified 
for the MFV case in Eqs. (39-37). Some loop corrections to the Higgs potential at large tan/3 
and their impact on tan (3 itself and the Higgs-fermion couplings have also been considered in an 
effective-field- theory framework in Ref. [58], which appeared during completion of this paper. 
Part of the results therein overlap with Section 3 of this paper. We disagree with [58] in some 
points (cf. Section 3), notably in that we find a tan /3-enhanced term in the relation of the DR 
and DCPR tan (3 parameters. We stress that, in general, only the former is numerically close to 
the tan P parameter extracted from S-physics observables. 

2. Large tan (3 phenomenology 

The prime application of our results is B — B mixing. We have identified a global U(l) symme- 
try of the hnQL Higgs-mediated FCNC transitions and the tree-level Higgs potential in the large- 
tan/? limit which suppresses the superficially leading contribution of Fig. 1. A systematic study 
of B—B mixing has required the analyses of four subleading contributions, which are governed 
by the small parameters mds/^6> 1/ tan/3, t;/MsusY and the loop factor l/(47r)^. These pa- 
rameters either provide a breaking of the C/(l) symmetry or allow for a contribution proportional 
to the C/(l) -conserving standard- model effective operator. Prior to this work, only corrections 
involving m,i jm,i) had been studied [5] (with the exception of Ref. [25]). The f /AfsusY correc- 
tions are found numerically small. The new loop contributions include all non-decoupling SUSY 
corrections to the quartic Higgs interactions A1-A7 and the contribution of neutral Higgs box di- 
agrams in the effective theory. In the complex MSSM the results for J^"^ comprising the neutral 
Higgs propagators become cumbersome. We have expressed in terms of sub-determinants 
of the neutral Higgs mass matrix. These expressions are easy to implement and clearly reveal 
the invariance of the Higgs-mediated amplitudes under rotations of the basis -sHd). The 
results for the Higgs sector are also used to refine the MSSM predictions for the S+ — > t^v 
and the Bg^a t^^l^" branching ratios. In this context we stress that loop corrections to the 
Higgs potential do not give rise to additional tan /3-enhanced contributions to the charged-Higgs- 
fermion couplings beyond those known before Ref. [25] appeared. Hence no modification of the 
charged-Higgs contributions to B^ t^v or B ^ Xg'y relative to Ref. [5] occurs. 

While the MSSM corrections to B—B mixing in the large tan (3 scenario could be dominated 
by the contribution of A5 and Ay, the size of this piece is limited by the experimental upper bound 
on B{Bs — > jjf^jjr). After performing an exhaustive analysis of this quantity, B{B — > Xg^), 
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— > T'^v) and the mass of the Ughtest neutral Higgs boson Mh, we find that the impact of 
the corrections to the Higgs potential on AMg is always weaker than that of the nis/ nif, correction 
identified in Ref. [5]. Assessing the total Higgs -mediated MSSM corrections to AM^ we find 
an upper limit of 7% of the SM contribution for > and Ma < 600 GeV. If fi is negative, 
the upper bound is around 20%. This is in contrast with Refs. [52] and Ref. [25], which claim 
large effects of the Higgs potential onB — B mixing. The corrections to B{Bq fJ-^fJ-^) from 
the Higgs potential are typically also small, but can reach 15% in some comers of the parameter 
space. In summary, the correlation between an enhancement of B{Bq — > n^ji') and a (moderate) 
depletion of AM^ found in Ref. [5] remains essentially intact. 

We finally note that our new contributions can alter the CP phase of the B — B mixing 
amplitude, while the previously known Higgs contribution proportional to mgrnhJ^^ has the 
same phase as the SM term (in MFV scenarios). While the maximal possible CP phase is clearly 
below the sensitivity of the current Tevatron experiments, it is an open question whether future 
B — B mixing experiments can help to unravel the CP structure of the MSSM Higgs potential. 



3. Heavy-quark relations and bag parameters 

We have transformed the NLO anomalous dimensions computed in Ref. [26] to an operator basis 
and a renormalization scheme typically used in lattice calculations. The according anomalous 
dimensions are needed to evaluate the 'bag' parameters, which parametrise the hadronic matrix 
elements, at the electroweak scale. We further employed a heavy-quark relation to improve the 
numerical prediction of the bag parameter Bf^^' entering the SUSY contributions to B — B 
mixing. The heavy-quark relation essentially determines Bf^^' in terms of the bag parameter 
BY^^, which is needed for the SM prediction [55]. We found Pf = -lBf^^'{mt) = -1.36 ± 
0.12. 
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A Notations and conventions 

To state our phase conventions for and M2 we quote the chargino mass matrix: 



( M2 

g V cos (3 
V V2 



V2 
1^ 



(98) 



with V — 246 GeV and the chargino mass term in the Lagrangian 



/rnass 



(99) 



For the case of a general flavour structure of the squark mass matrices we define the trilinear 
couplings Tu-^ and Td.^ (with flavour indices i, j) such that the squark mass matrices read 



Ml = 



V sin (3 



M~ 



T^- fj,*Y^ cot P 



V sin (3 



Mi 



UR 



I 



( 

V COS (3 



Ml 



Td-fJ,*YdteinP 



V cos P 



fJ-/.ritan/3 



(100) 



in the super-CKM basis, where the (DR-renormalized) Yukawa matrices are diagonal: Yq — 
diag (yg^ , y^j, y^j), q — u,d. The mass matrices in Eq. (100) correspond to the squark mass term 



d d 



(101) 



with $5 = {ul, cl, tL, ur, cr, IrY and $j = (d^, s^, 6i, (Ir, sr, hR^ . 

Our sign convention for the MSSM Yukawa couplings implies the following relations be- 
tween 2HDM quark masses and MSSM Yukawa couplings in the up sector: 



V2 

vsva.(3 



(102) 



In general, the analogous relations in the down sector involve complex phases associated with 
the tan /3-enhanced threshold corrections in Eq. (22). In particular, Ref. [24], which discusses 
the complex MSSM for the case without flavor mixing, relates the b Yukawa coupling to the b 
quark mass as 



Vb 



V2 

V cos /3 1 + 63 tan /3 ' 



(103) 
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rendering yi, complex for complex ?3 in Eq. (33). Our approach of matching the MSSM to an 
effective 2HDM permits different phase conventions, because the quark fields in the MSSM and 
the 2HDM can be chosen to differ by a phase factor. We can rephase the bp super-field of the 
MSSM in such a way that is real and positive and 1 + ?3 tan (3 in Eq. (103) is replaced by 
|1 + ?3tan/3|. The (physical) phase of 1 + 63 tan /3 will then, however, appear explicitly in 
the Higgs and higgsino couplings to bottom (s)quarks. Introducing 3x3 flavour mixing, the 
relation between Yd and m^, m^, rub is found from Eq. (22). Now the quark fields in the 2HDM 
differ from those in the MSSM by a complex rotation in flavour space and a particular choice 
for the phases of MSSM fields appears less obvious. In particular, one could render all y^^ real 
and positive by suitable rephasings of the right-handed superfields. Note that within MFV yh is 
still related to via Eq. (103) in good approximation without such rephasings. The analogous 
relation for the first two generations reads 

yd,s = ^T-, — 7 — ^, (104) 

V COS p 1 + eo tan p 

while in the lepton sector, we have 

Vi = ^ o 1 I — o with £ = e, T . (105) 

V cos p 1 + €£ tan p 

In (most of) the paper we express our results in terms of fermion masses (i.e. avoiding ydi,e) to 
achieve formulae which are independent of such phase conventions. Note that the phases of eo 
and ey are physical and no phase convention other than that of the CKM matrix matters for 
in Eqs. (29) and (30). While the phase convention of yg. enters the phases in Tg, it drops out from 
the MFV parameters ttg in Eq. (109). 

Finally, the quadratic squark soft-breaking terms are defined as follows: 

(Ml)ij = (VSkm rhlV^^KMh + ^-^^SM' + SijMl cos 2/3(1/2 - 2 sin^ Bw/S), 

cos^ 3 

(M|Ji,- = {m\)i^ + + 5iiM| cos2/3(-l/2 + sin^ ^iy/3), 



2 

{Ml^)ij = {rhDij + Y^5ij\yu,\'' + 2 (5i,M| cos 2/3 sin^ Ow/^, 



ip' cos'^ 3 



(106) 



where V^j^j^^ corresponds to the relative rotation of left-handed it-type and d-type quark fields 

performed when diagonalizing the Yukawa matrices. It differs from the actual CKM matrix, 
defined by the rotations that diagonalize the 2HDM mass matrices rather than the MSSM Yukawa 
couplings, by loop-suppressed (but tan /3-enhanced) corrections. In particular, within MFV, we 
have: 

^ + VSkm^, for (z, j) = («, 6), (c, 6), (t, d), (t, s), 

(107) 



VcKM.. = l l + ^3tan/3 



VniTM.. otherwise. 
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The relations Eq. (107) take a particularly compact form in the (exact) Wolfenstein parametriza- 
tion, where one has 



1 + eo tan (3 



A = A°, p = 77 = 77°. (108) 



1 + 63 tan (3 

Whenever we consider the case of MFV we write 
fuij = atVuAj, Td,j = abydAj, ^Qi^ = rn^q^ij, ml.. = ml5ij, ml.. = ml5ij. (109) 
The SU(2) relation between M~ and M~ then implies for the third generation: 

2 

777 

M| = M~2 + 777^ - - — ^ + M| cos 2/3(1 -sin^^vF). (HO) 

*^ ''^^ * |l + e3tan/3|2 z wj 

In the strict SU(2) limit (i.e., 7;/Msusy ^ 0) one has Mi = M? , but for small M? the term 

involving 777,^ can be relevant. Also FCNC W-uu loops vanish (for universal M~^) by the GIM 
mechanism up to the 777^ term in Eq. (110). 

Finally, it is convenient to define the so-called superflavour basis, obtained from a generic 
electroweak interaction eigenstate basis by rotating the supermultiplets Ql, ur and such that 
the quadratic squark soft-breaking terms are diagonal. We denote the corresponding entries by 
777q., 777^., and 777^.. For MgusY S> v, these are just the squark masses, and the computation of 
the effective couplings A, induced by heavy squark loops for arbitrary flavour and CP structure is 
greatly simplified. The Yukawa matrices and trilinear terms in this basis are simply written Y^^d 
and Tu,d, respectively. They are given in terms of Yu,d and T^^d as follows: 

^uSuV^j^^Vl Yj = UdYdV}, T^ = uXvSkmVI TJ = UdfdVl (111) 
where the matrices Vd, and Ud are defined such that 

diag(7fi|J = y,7f7|yi, diag(7n^J = C/„7f7^C/t, diag(7n^J = C/,7f7^C/J. (112) 

Assuming MFV, one is allowed to choose Vd = Uu = Ud = 1. 

Our conventions comply with the Les Houches accord [59]. In particular, our Y^^d and T^^d 
matrices correspond to a particular choice of the generic Yu,d and Tu,d matrices of Ref. [59]. 
Our conventions also agree with those of Ref. [60], except that the sign convention of our Yd in 
Eq. (103) is opposite. Besides, our equals — A„ and our Td equals Ad of Ref. [60], respectively. 



B Matching of the MSSM on a 2HDM 

The notation of Sect. 3 distinguishes between the coefficients Aj and Aj. The former quantities 
contain the results from the supersymmetric loop corrections to the quadrilinear Higgs couplings, 
whose tree-level values are given in Eq. (21). The latter coefficients also include the effect of the 
wave function and gauge coupling renormalization constants in Sect. B.l. 
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In the following we summarize the one-loop matching corrections for the quartic Higgs- 
coupling constants in the general MSSM. While the calculation of loop corrections to the Higgs 

sector of the MSSM has a long history, see for example [18-24], and a determined reader could 
extract part of the matching coefficients below from these works, the results collected in this 
appendix as a service to the reader are more complete than those in the literature, capturing the 
effects of the full set of mass, flavour-violation, and CP- violation parameters of the most general 
MSSM. 

The general results are quite lengthy hence we start with X2,5,7, where all renormalization 
constants are included, in the approximation of third generation dominance and degenerate soft- 
breaking parameters midQei — ^1,2 = ^1/2- The quartic coupling constant A2 is 
already present at tree level and as such depends on the renormalization scale /^o at one-loop. It 
reads: 



1 \r, ^, |4i 1|,,|4L, |4i 1|,,|4|,, |4 



4 167r2 

(6 \yt\ ^ - ^g^) \atyt\ ^ + Ivb] ^ + Ivr 



+ 



+ 5 



4 , log {M,) (13 - 3M^) nj^\ M, (5M^ + 14) + 1 

12 V /^o 



4(M^-1)3 



+ g g' 



2,./ 2 log (M^)(M^- 4) M2 (M^ + 5)M^ H/^ 



{M, - 1) 



(M,-l)2 
2 

2 log; "'2 



, 15, M^(17M^ + 158) + 5 , 

log (M^) (M^ ((141 - 43M^) - 12) + 4) 
12 (M^-l)3 



(113) 



where we have defined the mass ratio = iMlf^^/lf^l'^- Its renormalization- scale dependence 
is cancelled by the inclusion of electroweak corrections in the effective 2HDM. The other cou- 
pling constants important in the large tan /3 limit are 



A. 



1 J iJ^3al\y,\' + 3\yt\X + \yr\X) 



IGtt^ 



6m^ 



+ 



+ 2/52 + 3^/^^ (log + (log - 2) + 2) 



(114) 



(^,-1)=^ 
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and 

Ay- 



\ytt - j^fJ^g'^ {Gab IVbf - 3 \ytf at + 2 \yrf Gr)^ 

A* /o I |2 I |4 I I |2 /o I |4 I I |4 \\ 

(3|at| at\yt\ + |//| (Baftl^fel + a^j) 



6m* 



+ [{9' + 4,'^/ + 3,'^) ^ + 8 + 2,'^,^ + 3/) ^) 
(-M2 + 21og(M^)M^ + l) 



(115) 

X 



8(M^-1)3 

In the following subsections we quote the results for Aj in the general MSSM. The effective 

potential V in Eq. (12) must be used with Aj = Aj and the relation between Aj and Aj is given 
in Eq. (71); the renormalization constants needed in this relation are given in Sect. B.l. We 
decompose Ai_7 as 

\ino I \sferm 

A,=Af^ + . (116) 

The tree-level values Xf^^ are given in Eq. (21). A-'^" and Af'^'^™, given in Sect. B.2 and B.3, 
contain the contributions from higgsino and gaugino loops and from sfermion loops, respectively. 
Finally we also list the relevant loop functions in Sect. B.4. All these results are given in the 
superflavour basis including the most general soft-breaking terms. 



B.l Renormalization constants 

The renormalization of g^'^ is related only to the field renormalization of W and B, Zw,b — 
1 + SZw,B, if we decouple the sfermionic, higgsino and gaugino contributions: 5g' — —SZb/^ 
and 5g = —5Zw/2. The finite part of the one-loop wavefunction renormalization constants of 
the gauge bosons are 



6Zw — 



SZb^ 



9' 1 

iGTT^e 

9'" 1 

167r2 3 



4 log 



+ 81og^ + }_^ log^ + TVclo^ 



^2 



m 



Qi 



l4 



2Iog^ 



3 

E 

i=l 



log 



ml 



^1 



m 



2 + 9 log 2 



^Nc , Nc , m\ 



~ 2 



(117) 



where //q is the renormalization scale and the soft-breaking terms are written in the superflavour 
basis (see Appendix A). 
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The sfermionic contributions to the wavefunction renormaUzation constants of the Higgs 
bosons are 



327r2 



+ B'o{me,,mi.)Te.,T* 



ij (118) 

ij 

while the respective contributions of the gaugino and higgsino loops read: 

SZm = --^,1 {/Wi\M,\, + 3g'Wi\M,l |^|)) 

SZud = -jg^//* (/M^S^dMil, + 3/M*S;(|M2|, (119) 

B.2 Higgsino-gaugino contributions to A1-A7 

The situation of A5 = A5 is particularly simple: The matching correction only involves the box 
function and A™° can be written in a compact form: 

(Ml, M2, |/i|, + (120) 
///2m2l>o(Mi,Mi,|//|,|/x|) , 

if we use the loop functions defined in Sect. B.4. 
We find for A'"^" = A^, . . . Ar, A^, Aij^°: 

A^°° ^g^ (a, + a2£'2 (M,, M^, H |//|) + a^D^ (M2, M,, |//|)) + 

g^g'" (< + 4^2 (Ml, M2, + a'^D^ (Mi, M2, + (121) 

/ (< + a'^D2 (Ml, Ml, I/.I) + a'ii), (Mi, Mi, |/.|, |/.|)) , 
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where the coefficients ag, . . . a'l depend on the index i labeUng Aj in Eq. (116), which we suppress 
throughout this appendix. The coefficients 03 are given in Table 4, while 

0^ , a'^= I 1 , a'l = \ 0^ for { A4 (122) 




A 


02 


04 


a'2 


04 


4 




Ai 




5 
2 


i(MiM* + M*M2) 


1 


ilMip 


1 

2 


A2 




5 
2 


i(MiM* + M*M2) 


1 


ilMip 


1 
2 


A3 


3|/x|2 + f|M2p 


1 

2 


2|/x|2 + l(MiM* +Mi*M2) 


1 


|^P + i|Mi|2 


1 

2 


A4 


-3|/ip-2|M2p 
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-2 
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fi (Ml + M2) 





/iMi 





At 


3/XM2 





II (Ml + M2) 





^Mi 






Table 4: Coefficients entering A^-A^"", A|?°, and X\^° in Eq. (121). 



B.3 Sfermion contributions to A1-A7 

The sfermion contribution to Ai_7 are products of loop functions and flavour dependent coeffi- 
cients if we sum over the generation index of the internal sfermions. For A5 our results then take 
the simple form 

^sferm ^f^vkijj^ [rfie^, m^j , m^, rfii^) + d'^^D^ {rhdi,rhdj,fhQ^,rhQ^) + 
d{ Do [niQ. , mq. , , m„ J , 

where the slepton contribution is contained in d*/'^' listed in Table 5 and d'^^l comprises the 





dT' 


Ai 




A2 


-\fl\^Ye,.Ye,-Y* Y* 
\r~\ ^kj efci eij 


A3 


I//|2t^ V I T* V* -I- T* V* \ 


A4 


\a?%,T: Y,, Y; 


As 


-fX^Te,.Te,-Y* Y*. 


Ae 




At 


^^\^^\ TeijYe^.Y*^y*^, 



Table 5: Slepton contribution to Af >-"^-Af ''"^ in Eqs. (123), (125), and (128). 
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4^'' 




Ai 
A2 

A3 
A4 
A5 
Ae 
At 


~-^-'-dki-'-dij-'-dkjdu 

3M'Td,Yd,Yiyi 


— '^-LUik-LUjl-Luil-^Ujk 

3 Tu.iT*.^Yuf^Y*.^ 

— 3/x Tu.j^Tu^jFJ.^F^*^ 



Table 6: Do squark contribution to Xf^"^-Xf^^ in Eqs. (123), (126), and (128). 



squark contribution (Table 6). Only A4 receives a contribution from 0^4 



df' = -3(r,,T4 - \fi\'Yd,Y:jiV* - \fi\'Y^,,Yn for A4, (124) 



while rfy^' = for Aj with i ^ A. The contributions to the matching coefficients depend on the 
Yukawa couplings Y^^u^d of the charged leptons, the up-type quarks, and the down-type quarks 
as well as on the trilinear soft breaking terms Tg^^ d, defined in the superflavour basis (see Ap- 
pendix A). 

We write Afl™ = A^L4 + Ai'L4, and find for the slepton contribution 



{b4Sij + hY^eiAj + beYeeiYeeji) Bq {fki^^fhi.) + 

{c,\i,\'y,y:j,, + c,t,,t:,m + c7i/^i%,nLne.,+ 

c^Te^iTe^Yeek,) Cq {rhe,, rhi^ , rhi^) , 
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while the squark contribution reads 



^^1-4 = (^7% + hYddJij + hYdd^Ydd,,) Bo {rhdi.rhd^) + 
hioYdu,yudj,Bo {rhdi,7huj) + 
(^11% + buYuu.Aj + bi3Yuu,yuu,,) Bo [fhu.^rhuj) + 
(6i4% + bi^Ydd^Aj + bieYuu.Aj + &i7'fdd,,'?dd,i+ 

bisYuuijYuuji + bigYddiyuuji) Bq (niQ^, rriQj) + 

(cM'Yd,y:j,j + c,oTd,raJi3 + ciil/^l 



{c,M^Yd,y*j,k + ci4|//r%i7,/d<i., + ci5r,,,r*/,-,+ 
(^-ci8Yd.y*.YduM'^ - ci8r„.,r;^,T„d,j/i|V 

{cM'Y^,,Y:J,j + cM'Y^,,YlY^u,, + C2iTu,XJij+ 

C22TuikT*,Yuu^. + C23Tu^^T*.Ydd^^^ Co {mQ.,mQ.,mui^) + 

{c2M'Y^,y:jjk + cM^Yu^^Yiyuu,, + c26t„,t;/,,+ 

C27T^,TlYuuj^) Co {mQ„rhuj,rhu,) + 
d'l^^Do {rhe^ , TJiej , fhi^ ,'ml^)+ d'f^Do {rhd^ , rhd^ , rhq^ , mg, ) + 

d'i^^Do {mQ,,mQ.,muk,mui) + d'f^Do {rhdi,rhQ.,rhQ^,rhu,) . 



Here we introduced shorthand notations for the products of two Yukawa coupling matrices: 



Y^Vij — Yxii^vij ' Y^Vij — Yxii^yij ' (127) 



where x,y — e,u,d and we sum over the internal index I. The coefficients 6„, c„, and d^^^^ for 
each Ai , . . . A4 are given in Tables 5, 6, 8, and 9. 
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\ (-/ - 3^^) 
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J(3S^-/) 
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-9' 
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Table 7: Coefficients of Af ™ and Af ™ in Eq. (128). 



We finally give the slepton and squark contributions to X^f. 



[c!^HTe^Y*..5jk + C^Te-.y^yee^J C'o {rfle^, TTlij , rfii^) + 

(c'jnTd.Y^Jjk + c'sfiTd-YlYdd^^^ Co {md^,mQ.,mQ^) + 
(c'QiJ,T^.^Y*Jij + c[ofiTu^^Y*^Yuu,,^ Co {mQ.,mQ^,mu,) + 

(c[^fJ^TuiY*..djk + c'i2fJ'TuiyiYuujk^ Co {mQ^,rhuj,rhuk) + 

d^'^^Do {rfie, , rfiej , m/^ , m/,) + 4'^'''i:)o (?Tt<i, , rfidj , mg^ , mg,) + 
4^'''i:)o (mg, , niQ. , rfiu^ , ) 



(128) 



where the coefficients c'^ and d^'^' are given in Tables 5, 6, and 7. 
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Table 8: Slepton loop contributions to Af 



sferm 



\ sferm 
A4 



in Eq. (125). 



B.4 Loop Functions 



In the UV-divergent loop functions we set e = {A — D)/2. The loop functions are defined as 



-,Aq (mi) -^e 



(47r) 



;Cq {nil, 1712,1713) 



(47r)- 



;Do (mi, m2, ma, 7714) 
-ly (mi,m2) 
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47r 
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Att 



-IE 



-lE 



-lE 



(27r)^ g2 



mt 



(27r)^ g2 
d^q 



m\ q^ 



{27r)^ — mf q^ — m\ q^ — m\ 



d^'q 



(2 



-R}^ q^ — m\ q^ — rn^ (f — vn\ (f — m| 



d 
d]^ 



fe2=0 



{27r)D {{q - ky - ml) {q^ - ml 



(129) 
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Table 9: Squark loop contributions to Af^'^™ . . . A| ™ in Eq. 



(126). 
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These functions read: 



Aq (mi) ^^+ml + ml log ( 



Bo (mi, 777,2) = - + 1 + 



Bq (mi, ma) 



m^log (^)+milog (^) 
m| — m2 

4 4in2 2i /'T 

mi - m2 + 2mi m2 log 



2 (mf — m|)^ 



m^mpogf^j +m^mpog(^) +m?m|log'^ 



Co (mi, m2, ms) = . o oVT^ 2^ ~ ^^^^^ 

[mf — ma) (mj — m^) (ma — m^) 



{mj,m2,m3,m|} 
+cyclic permutations 



Do (mi, ma, ma, m4) = _ (n-A (n-^^ (h - A (h - (r-rl' 



{a — b) {a — c) {a — d) {b — c) {b — d) (c — d) 



W-(mi,m2) = ---21og f4 



log -I 



(2m2 — Qm\m\) m\ — &m\rn\ + ma 



m: 

^1/ (mi— ma)^ (mf — m^)^ 



-Da(mi, ma, ms, m4) =Co(ma, ma, m4) + mlDo{mi, ma, ma, m4) 
£)4(mi, ma, ma, m4) =5o(m3, m4) + (mi + ma) Co(ma, ma, 1714)+ (131) 

mi Do (mi, 1712, ms, rn4) 

A further loop function, H2, is defined in Eq. (34). 

C Renormalization group and bag parameters 

The standard-model contribution to 5— S mixing involves the operator Qy'"'" = {bLliJ.(lL)ibLl^(lL) 
of Eq. (6). The main new supersymmetric contribution io B — B mixing presented in this paper 
comes with the four-quark operator — (bRqL)(bRqL) with q = d or q = s, see Eq. (2). 
qSll jnixes under renormalization with 

Qf" = (M(^«?l) (132) 

where 7, j are colour indices. The operators Qf^^ and Qf^^ are widely studied in the context of 
the width difference AF among the two mass eigenstates in the B—B mixing system and the CP 
asymmetry Ofg in flavour-specific decays [61,62]. 

Yet the next-to-leading-order (NLO) anomalous dimensions have been calculated for an 
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equivalent operator basis in Ref. [26]. These operators, 

Q2^ - ibRqL)ibLqR) , 

Qf"^ = {hnqL^b^L) , 

Qf'^ = -{bRa^,qL)(bRa^''qL), 

Q7^^{bRl,qR){bRl''qR). 

QT^ = {pLqRKbLqR) , 

Qf^ = -{hcT,,qR)iha'^''qR), 

are split into five sectors (VLL, LR, SLL, VRR, SRR) which separately mix under renormaliza- 
tion - note that we define a^^ = | [oV-7i/]- The anomalous dimensions of the VRR and SRR 
sectors are the same as those of the VLL and SLL sectors, respectively. To define the renormal- 
ization scheme for the NLO we first note that we use the MS scheme with anticommuting 75 as 
in [26]. Then we must specify the definition of the evanescent operators which enter the NLO 
results as counterterms. In particular for the SLL sector the evanescent operators of Ref. [26] 
read: 



SLL _ ^ui „ „3\ (U „IJ-i'„i\ «/lSLL 1/^SLL 

(134) 



SLL 

2 : 



Ef^ = (&'r7m7.7.7.?1) (J^^rI'^YYiWl) + (-64 + 96e)gf ^ + (-16 + %e)Q 
Ef^ = (^^7M7.7p7agi) {yRrrrrqi) - 64Qf ^ + (-16 + 166)gf L , 

where we use e = {A — D)/2. 
The operator basis 

Q2 — Q2 -I 

gSLL ^ gSLL ^ 

QfL = QfL =(6^gi)(6X), ^^^^^ 

gVRR ^ gVRR ^ 



which we adopt in this work agrees with the one of Eq. (133) except for the SLL sector and the 
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SRR sector. The evanescent operators are defined as in Refs. [55,61]: 

={Al,l^<fLWRl''iyi) + 8(1 - e)gf ^ , 
' HhRl.luq'^MyRYY^L) + 8(1 - 6)gf ^ . 

The hadronic matrix elements in this basis are parametrized in terms of 'bag' parameters 
^VLL^ ^SLL/^ ^SLL/ defined as 

{B,\Qf\^)\B,) = -Am^^ /I ^'(/i), (137) 

{B,\Qf\„)\B,) = ^MljlBf^^\„). 

Here fj, is the renormalization scale at which the matrix element is computed and /^^ is the Bq 
meson decay constant. While fs^ exceeds fs^ by 10-30%, no non-perturbative calculation finds 
any dependence of a bag parameter on the flavour of the light valence quark. In the vacuum 
insertion approximation the bag parameters equal BY^^{ij) — 1 and B^^'{ij) — Bf^^'di) — 
^Bq/V^bilA + ^qilj)]^- Lattice computations determine the matrix elements at a low scale 
around 1 GeV and results are quoted for /i = mf,(mb). In order to use the lattice results in our 
calculation we need the renormalization group (RG) evolution of the bag parameters to the high 
scale jih which is set by the masses of the Higgs bosons exchanged in our B—B mixing diagrams. 
The matrix elements computed on a finite lattice are converted to continuum QCD by a match- 
ing calculation. This lattice-continuum matching is only meaningful beyond the leading order 
of perturbative QCD. Thus the dependence of the bag parameters on the chosen (continuum) 
renormalization scheme must be addressed: The NLO anomalous dimension matrices entering 
the RG evolution must be defined in the same renormalization scheme as the bag parameters, 
so that the scheme dependence properly cancels from physical observables. The NLO anoma- 
lous dimensions have been calculated for Qj^^ in Ref. [63]. As said previously, in the case of 
(Qi^^, Q^^) the NLO anomalous dimensions have been calculated for the equivalent operator 
basis (Qi^^i Q'¥'^) with the evanescent operators of Eq. (134) [26]. 

The purpose of this^ection is twofold: First, we present the transformation of the results of 
Ref. [26] to the (Qi^^, Q^^) basis and the scheme corresponding to the evanescent operators of 
Eq. (136), for which lattice groups quote their results. These formulae are useful beyond the need 
to evolve the bag parameters given at // = up to /x = /x/i : In particular lattice groups need to 
evolve Bf^^'di) and Bf^^'^i) from a scale around 1 GeV up to = mj. Second, we exploit a 
heavy-quark relation among the bag factors in Eq. (137) to sharpen the numerical prediction for 
Bf^^'(Hh) entering the SUSY contribution to B — B mixing. While constraints from the heavy- 
quark limit of QCD have been used to improve the predictions for AF and Ofs [55, 61, 62], they 
had escaped attention in studies of new physics contributions to B physics observables so far. 
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C.l NLO scheme transformation formulae 

We decompose the anomalous dimension matrix in the usual way as 



7 = 



M/f) (0) 

47r ' 



An 



7^^) + 0{al] 



(138) 



The NLO correction 7^^) has been computed for the basis {Qf'^, Qf'^) in Ref. [26]. In four 
dimensions it is related to the basis (135) by a simple Fierz identity: 



Q 



D=4 



R 



gfL 



RQ, 



(139) 



where R is given in Eq. (141) below. 

Yet in D dimensions our change of basis involves a rotation of the operator basis - including 

the evanescent operators E = {Ef^^, E|^^)^ 
follow Ref. [64] and write the rotation as ^ 



{Ef^'", E2^^y - and a change of the renormalization scheme. We 



Q ^ R(^Q + WEy E^M[eUQ + 



with 

R 



I 
1 1 

"2 8 



w 




8 



U 



1 _ : 
4 'i 
8 



1 + 



M 



E 



)■ 



-8 
-4 1 



(140) 



(141) 



The information on the definition of the evanescent operators in Eqs. (134) and (136) is con- 
tained in the matrices U and M. Now Eq. (140) corresponds to a finite renormalization with 
renormalization constants [64] 



^(1,0) 



R 



WZ 



(1,0) 
EQ 



^ + wz^^ - t^^m ) u 



1 



R- 



(142) 



While the one-loop anomalous dimension matrix is just rotated, the two-loop anomalous dimen- 
sion matrix undergoes an additional scheme transformation: 



7W = ^7(°)^-\ 
7«=^7«^-i- 



^QQ ,1 



2^(0)^(1,0) 



(143) 



QQ ' 



with the one-loop operator renormalization constants 



7(1,1) _ 



-8 -8 



Mi,i) 

'EE 



16 
3 



n 

3 



^(1,0) 
^EQ 



37 



29 

48 



12 



(144) 



'Two more evanescent operators (called i?f and Ef^^ in Ref. [26]) must be specified to fully define the 
scheme of the calculated 'y'-^-'. This information enters the matrix U in Eq. (141). We choose to add E^^^ and 
£^4 also to the evanescent operators of Eq. ( 1 34), so that we can in practice work with the change of basis defined 
in Eqs. (140) and (141). 
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We can now calculate the new two-loop anomalous dimension matrix 7 from the NLO 
anomalous dimension matrix 7 of Ref. [26], 



We obtain 



-10 - 



-40 



7 



(0) 



V 



1 

6 

34 

y 



28 4 \ 

y 3 

16 32 



7 [26] 



1459 74 _35 _ 1 \ 



6332 584 2065 



394 



27 y 



(145) 



7 



(1) 



V 



260 88 44 8 \ 
"~3~^27'^ ~y^27'' 



3 27-^ 27 



(146) 



Here / denotes the number of active flavours and 7^'^^ coincides with the result in [61]. As a 
check we have calculated the result of Eq. (146) also in a different way: It is possible to define 

evanescent operators such that the Fierz identity holds for the one-loop matrix elements. This 
choice fixes the definitions of both Ef^^ and i^f^^ in Eq. (136) and of the evanescent operators 
on the {Qf^^, Ql^^) basis. (One of the latter operators equals e times a physical operator. Its 
impact is equivalent to a finite multiplicative renormalization of Qf^^.) In this approach one can 
simply rotate 7^^^ in the same way as 7^°^ in Eq. (143). Finally the result is transformed to the 
scheme of Ref. [26] using the scheme transformation formula of Ref. [65]. 

Next we calculate the matrices governing the RG evolution in the [Qf^^, Qf^^) basis. The 
bag factors at the scale /x/j are obtained from those at the low scale = 0(mb) via 



= U{iib,fJ'hf 



(147) 



In the spirit of [27] we write the evolution matrix as 



(Hb)^ ' 47r ^asii^b) 



where U^^^ is the LO evolution matrix and the NLO correction reads 



(148) 



(149) 



The 2x2 matrix Jf is calculated from the anomalous dimension matrix 7 [66]. We only need 
J5, since we run with 5 active flavours to the scale For applications in kaon physics one also 
involves J4 and J3. We quote all three matrices here, so that the formulae of Ref. [27] can be 
easily extended to the {Qf'^, Qf'^) basis: 



J. 



1.474 0.707 
0.306 -5.350 



J4 



0.964 1.452 
0.375 -4.982 



0.652 2.597 
0.421 -4.804 



(150) 
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We quote handy formulae for the five-flavour evolution matrix, similarly to Ref. [27]: 

.7(0) _ f 0.9831 -0.2577 \ 0.6315 , 0.0169 0.2577 \ o.7i84 



f=5{V) = 1^-0.0644 0.0169 J ^ ^ \ 0.0644 0.9831 J ^ 
The NLO correction reads: 

/ 1.4040- 1.3707^ -0.3680 - 2.073177 \ 0.6315 

""f-^^''^ \^ 0.6454 + 0.0898 ?7 -0.1692 + 0.1358 J ' 

f 0.0704 - 0.1037?7 1.0746 + 1.3665 r/ \ o.rm 
^ V -0.3395 - 0.3958 -5.1807 + 5.2141 )^ ' 

In our numerical analysis we drop the terms which are linear in rj in the two matrices in Eq. (152), 
because they are scheme-dependent. The scheme dependence of these terms cancels with that of 
the NLO QCD corrections to the 5—5 mixing diagrams with SUSY Higgs exchange. Yet these 
QCD corrections are unknown. 

C.2 Hadronic matrix elements and heavy-quark relations 

The three bag factors B\^^, Bf^^'dj.^), and Bf^^'(ij,i,) obey a heavy quark relation [62]: 

^^"'(a^.) = ^«.(A.,)Sr + l^MBf'^''' + O . (153) 

Here q;i(//) and q;2(a*) comprise NLO QCD corrections [55,61]: 

47r \ mi, J Att \ nib 6 J 

These values are specific to the definition of the evanescent operators as in Eq. (136). As men- 
tioned in Sect. C.l, this definition allows to maintain the validity of Fierz identities at the loop 
level. Such a definition is preferred, if the bag factors are meant to parametrize the deviation 
of matrix elements from the vacuum insertion approximation (VIA), because the calculation of 
matrix elements in VIA approximation involves a Fierz transformation. In particular the choice 
in Eq. (136) is crucial for Eq. (153) to hold in the limit of a large number Nc of colours [61]. 

The bag factor BY^^ is very well studied in lattice QCD, so that it is worthwile to study the 
constraint on the other bag factors when Eq. (153) is combined with lattice results for B^^^. 
Indeed, one can use Eq. (153) to pinpoint the ratio 

J.,,,: , = 0.93 + 0.23 J.,,,) { + 0.23 + 0.05 — - (155) 

quite precisely, even if Bf^^' is only poorly known, because its coefficient in Eq. (155) is small. 
The last term in Eq. (155) quantifies the Agco/^^b corrections, see [55] for details. The lattice 
results of [67] have been combined in Ref. [55] to 

BY^^irrib) = 0.85 + 0.06 and Bf-^'irub) = 1.41 + 0.12. (156) 
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Inserting these values into Eq. (155) yields 



1.57 ±0.08, (157) 



which is consistent with the direct determination 

Bf^^'irub) = 1.34 ±0.12 (158) 

from the lattice [67]. 

We are now in the position to accurately predict the bag factors at the high scale Hh- Choosing 
= mt(mt) = 164 GeV, as(Mz) = 0.1189 and mb{rnb) = 4.2 GeV and using Eqs. (157) and 
(156) we find 

Bf^^\mt) = 1.62 5fL'(m5) +0.01 5fL'(mb) 

= (2.54 ± 0.13) Bi^^^(mfe) + 0.01 
Bf^'^\mt) = 1.29Bf^^\mb) + 0.5ABf^^'{mb) 

= (2.03 ± 0.10) ^^^^(mb) + 0.77 + 0.07 (159) 

Here we have omitted the scheme-dependent terms proportional to 77 in Eq. (152). The small 

(2, 1) element of uflr^ in Eq. (151) ensures that Bf^^'{mi,) is inessential for Bf^^'(mt). One 
realises from Eq. (159) that the uncertainty of the high-scale bag factors stems almost completely 
from the error of the lattice result for BY^^{mb). 

Switching finally to the P/s defined in Eq. (8) we get 

pSLL ^ _^Bf^i^'{mt) = -(1.59 ± 0.08) 5y^^(mb) -0.01 = -1.36 + 0.12 

8 

pVLL ^ B^^^{mt) = 0.83 B^^^irrib) = 0.71 + 0.05. (160) 

In the last line the full NLO result of [63] has been used. We don't need P^^^ = Bf^^ ' (mt) /8 for 
our analysis. Parity ensures that and the chirality-flipped operator defined in Eq. (7) 
have the same matrix element, i.e. P^^rr _ pSLL 

Finally we compute P^^ using the formulae of Ref. [27] with the bag factors of Becirevic et 
al. [67]. This time the conversion between the bases of Ref. [27] and Ref. [67] is straightforward, 
since the renormalization scheme used in Refs. [26, 27] respects the Fierz symmetry and lattice 
results are already quoted for this scheme. The result is 

P^^ = 3.2 + 0.2. (161) 



The number in Eq. (161) is significantly larger than Pg'"^ = 2.46 quoted in Ref. [27], because our 
value for vrib is smaller and the lattice bag factors are larger than one. The error in Eq. (161) does 
not include the systematic error from the quenching approximation. 
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D Trilinear Higgs couplings 

The trilinear terms of the effective Lagrangian at tan (3 — oo introduced in Sect. 2.3 read 



Again, the first two lines respect the f/(l) symmetry introduced in Sect. 2.3, while the last two 
lines break it, and the breaking is proportional to loop-induced couplings. Finally, the quartic 
Lagrangian is obtained from the quartic terms in Eq. (12) by substituting if„ — > (/t+, -^<Pu) ^^'^ 
Hfi ih^d-:~^d )- ^^^^ there, only A5, \q and Ay break the symmetry. 
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